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Abstract 

This  study  performed  an  investigation  on  determining 
the  appropriate  number  and  locations  of  continental  United 
States  aerial  ports.  To  accomplish  this  a  linear 
programming  formulation  was  adapted  with  the  optimizing 
function  based  on  trading  off  the  cost  of  shipping  cargo 
against  port  operating  costs.  Cargo  would  travel  from  CONUS 
origin,  through  aerial  port  of  embarkation  (APOE) ,  to  aerial 
port  of  debarkation  (APOD)  at  minimum  cost  to  the  DOD.  The 
need  for  the  study  was  precipitated  by  continued  reductions 
in  the  military  budget,  consolidation  of  defense  depots,  and 
the  reduction  in  the  number  of  personnel  stationed  overseas. 

Cargo  movement  data  was  extracted  from  the 
Transportation  Reporting  and  Inquiry  System  database  for 
fiscal  year  1996.  This  information  was  then  used  as 
deterministic  demand  at  the  APODs  from  particular 
origination  cities i  The  demand  had  to  be  exactly  met  in  the 
formulation.  Applying  the  linear  program  resulted  in  the 
recommendation  to  operate  only  three  aerial  ports.  They  are 
Travis  AFB,  CA,  Dover  AFB,  DE,  and  McGuire  AFB,  NJ  saving 
over  11  million  dollars  a  year. 
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AERIAL  PORT  LOCATION  STUDY 


I .  Introduction 

At  a  time  when  the  United  States  Armed  Forces  is 
reducing  its  numbers  in  response  to  a  changing  world 
situation  and  budgeting  constraints,  there  is  a  need  to 
readdress  the  structure  of  the  worldwide  cargo  distribution 
system.  The  military  is  now  emphasizing  continental  United 
States  (CONUS)  based  forces  with  the  ability  to  redeploy  on 
short  notice.  Of  course,  with  fewer  personnel  in  forward 
areas,  the  total  supply  requirement  will  decrease.  This 
reduction  will  precipitate  a  restructuring  of  the 
transportation  system. 

A  major  portion  of  this  system  consists  of  aerial  ports 
of  embarkation  located  throughout  the  United  States.  When 
the  aerial  port  distribution  network  was  originally 
designed,  aircraft  were  of  limited  range  and  capability  as 
compared  to  today's  standards.  It  made  sense  in  the  post 
World  War  II  environment  to  locate  most  of  the  aerial  port 
facilities  along  the  east  and  west  coasts.  But  today  with 
long-range  transport  capability  of  United  States  Air  Force 
(USAF)  and  commercial  aircraft,  a  new  look  at  the  basing 
structure  is  warranted. 
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In  1989  a  study  along  the  lines  of  this  thesis  was 
conducted.  It  was  called  the  "Optimal  Airlift  Distribution 
Study  Proposal"  (OADS)  completed  by  Greg  Holevar,  currently 
at  headquarters  Air  Force  Material  Command  in  the 
transportation  directorate.  The  OADS  study  had  some 
interesting  results.  At  that  time,  the  current  locations 
considered  were  the  eight  major  aerial  ports  plus  two  inland 
bases.  They  were  Charleston  AFB,  Charleston,  SC;  Dover  AFB, 
Dover,  DE;  McChord  AFB,  Tacoma,  WA;  Norfolk  NAS,  Norfolk, 

VA;  Norton  AFB,  San  Bernadino,  CA;  Tinker  AFB,  Oklahoma 
City,  OK;  Travis  AFB,  Fairfield,  CA;  McGuire  AFB, 
Wrightstown,  NJ;  Kelly  AFB,  San  Antonio,  TX;  and  Hill  AFB, 
Ogden,  UT.  After  consideration  of  cargo  origin, 
destination,  and  aerial  port  effectiveness  in  handling  the 
calendar  year  (CY)  1988  cargo,  the  study  recommended  a  few 
changes.  The  optimal  locations  of  ports  were  determined  to 
be  Charleston,  Dover,  Norfolk,  Tinker,  and  Travis.  It  was 
recommended  that  McGuire  and  McChord  downsize  for  a  wartime 
role  and  that  Norton  close  its  doors  (1)  . 
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Specific  Problem 


Since  1989,  at  the  completion  of  the  OADS  and  its 
implementation  beginning  in  1990,  there  have  continued  to  be 
drastic  reductions  in  personnel  and  base  infrastructure. 
These  are  the  results  of  continued  tightening  of  the 
Department  of  Defense  (DOD)  budget  and  Base  Realignment  and 
Closure  Commission  (BRAC)  impacts.  Bases  overseas  and  in 
the  CONUS  are  still  being  shut  down  or  realigned  to  include 
the  depots  and  air  logistics  centers  (ALCs) ,  from  which  most 
of  the  cargo  originates  that  enters  into  the  Defense 
Transportation  System  (DTS) .  As  a  result  of  BRAC  95, 
several  depots  are  slated  for  closure  by  the  year  2001.  The 
current  list  of  Defense  Distribution  Depots  or  Facilities 
contain  24  sites.  Those  depots  scheduled  to  close  in  1997 
include  Letterkenny,  PA,  Ogden,  UT,  Columbus,  OH,  and 
Memphis,  TN.  The  depot  at  McClellan  AFB,  CA,  and  the  San 
Antonio  ALC  (SA-ALC)  at  Kelly  AFB,  TX,  are  scheduled  for 
closure  in  2001  (2) .  As  the  independent  service  depots 

continue  to  merge  into  the  Defense  Logistics  Agency  (DLA) , 
the  number  of  cargo  origination  points  will  constrict  even 
further  necessitating  another  look  at  port  location. 

Today' s  environment  is  a  reflection  of  a  similar 
situation  in  1989  and  1990.  As  taken  from  Annex  Alfa  to  MAC 
(Military  Airlift  Command)  Programming  Plan  90-16,  this 
quote  sums  up  the  situation: 
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Current  and  future  [DOD]  budget 
constraints  require  a  closer  look  at  the 
current  way  of  doing  business.  The 
existing  and  future  constraints  on 
Second  Destination  Transportation 
funding. .. [and]  MAC'S  task  of  meeting 
the  time  standards  set  by  the  Uniform 
Material  Movements  and  Issue  Priority 
System  (UMMIPS)  is  becoming  an 
increasing  difficult  and  expensive 
challenge.  (3) 

According  to  Mr.  Stef fey  (4),  at  HQ  AMC  Cargo 
Management,  approximately  50  percent  of  the  cargo  handled  by 
AMC  is  currently  not  meeting  the  Uniform  Material  Movement 
and  Issue  Priority  System  (UMMIPS)  requirements.  As  a 
result  of  this  and  other  tumultuous  changes  occurring  within 
the  DOD,  Air  Force  Materiel  Command  (AFMC)  and  Air  Mobility 
Command  (AMC)  determined  that  it  was  again  time  to  address 
the  DTS  structure  and  more  specifically  the  basing  strategy 
of  the  CONUS  aerial  port  system.  The  results  may  be  the 
more  efficient  and  effective  handling  of  the  millions  of 
shipments  that  cross  AMC's  path  annually. 


Research  Scope 

This  thesis  will  analyze  the  cargo  data  provided  by  the 
Transportation  Reporting  and  Inquiry  System  (TRAIS)  (5)  to 
determine  the  origin  and  aerial  port  of  debarkation  of  the 
majority  of  the  cargo  that  transit  the  DTS  system  for 
overseas  delivery  via  AMC's  aerial  ports.  The  purpose  of 
this  effort  will  be  to  locate  a  best  set  of  CONUS  aerial 
port  facilities.  The  set  of  alternatives  include  six  major 
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aerial  ports  in  existence  today  along  with  three  interior 
possibilities . 

Using  this  synthesis  of  TRAIS  information,  a  linear 
programming  algorithm  will  be  developed  to  enumerate  and 
locate  a  possible  set  of  aerial  port  facilities  to  handle 
the  Second  Destination  cargo  requirement. 

Model  formulation  will  be  based  on  a  transshipment 
facility  location  problem  which  locates  intermediate  service 
facilities  to  minimize  some  objective  cost  function,  or  in 
this  case,  the  total  cost  of  shipping  cargo  through  the 
aerial  port  system  en  route  to  the  final  destination 
considering  operating  costs  of  the  aerial  ports. 

Research  Questions 

There  are  two  questions  which  this  thesis  will  answer. 

1.  Based  on  a  modified  distribution  system  location 
problem,  what  are  the  optimal  locations  and  how  many 
aerial  port  facilities  are  needed? 

2.  How  does  the  new  structure  compare,  on  the  basis  of 
cost,  with  the  current  APOE  structure? 

Chapter  III  will  discuss  the  methodology  used  to  answer  the 
above  questions. 

Summary  and  Overview 

This  chapter  has  given  a  brief  introduction  and 
background  to  the  subject  of  aerial  port  facility  location. 
It  has  also  addressed  the  specific  problem,  research  scope. 
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and  the  questions  that  will  be  answered  as  a  result  of  this 
effort.  Chapter  II  will  review  the  literature  expounding  on 
location  problems  in  general  on  down  to  the  focused  task  at 
hand.  The  methodology  will  be  the  thrust  of  Chapter  III, 
which  will  yield  the  results  in  Chapter  IV  and  conclusions 
in  Chapter  V. 
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II.  Literature  Review 


Facility  location  problems  have  been  the  subject  of 
research  for  quite  some  time.  Advances  in  computing  power 
in  the  last  10  to  15  years  have  caused  a  resurgence  in  the 
genre.  For  excellent  discussions  of  the  basic  facility 
location  problems  and  models,  three  sources  come  to  mind. 
The  first  is  a  textbook  on  discrete  location  theory  by 
Mirchandani  and  Francis  (6)  and  the  second  is  a  text  on 
facility  layout  and  location  analysis  by  Francis,  McGinnis, 
and  White  (7).  Another  good  treatment  of  the  subject  was 
authored  by  Love,  Morris,  and  Wesolowsky  (8). 

Mirchandani  and  Francis  narrow  their  focus  to  discrete 
location  decisions  as  opposed  to  continuous  location 
decisions : 

The  major  reasons  are  that  in  most  cases 
decision-makers  consider  a  discrete 
representation  to  be  a  more  realistic 
and  a  more  accurate  portrayal  of  the 
problem  at  hand,  and  that  continuous 
formulations  appear  to  be  relatively 
difficult  to  solve.  (6) 

They  present  formulations  of  and  solution  methods  for  the 
basic  models  and  their  variations.  These  would  include  the 
p-median,  the  p-center,  the  uncapacitated  facility  location, 
and  the  quadratic  assignment  problems  (6). 

Francis  et  al.  discuss  both  planar  single  and  planar 
multifacility  location  problems.  They  also  describe 
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transportation  network  setups,  to  include  the  tree,  median, 
center,  covering,  and  warehouse  location  problems  (7). 

Facility  location  models'  main  purpose  is  to 
quantitatively  evaluate  the  alternatives  of  siting 
facilities,  be  they  warehouses,  plants,  etc.,  to  minimize 
the  cost  or  some  other  objective.  Models  are  mathematical 
optimization  techniques  that  can  be  used  to  determine 
whether  or  not  a  facility  should  be  opened  or  closed,  and 
where  they  should  be  placed  (9). 

The  numerous  assumptions  made  that  simplify  any 
particular  application,  ultimately  determine  the  solution 
generated  by  the  model.  One  assumption  in  many  location 
models  is  that  the  demand  to  be  satisfied  by  particular 
facilities  are  fixed  and  known  (9) .  Estimates  of  capacities 
and  costs  are  also  used.  The  costs  may  be  divided  into 
transportation,  fixed,  or  operating,  with  the  assumption  of 
linearity  for  the  transportation  portion.  Therefore,  the 
accuracy  and  quality  of  results  are  heavily  dependent  on  the 
realism  associated  with  these  assumptions. 

Another  point  to  be  made  about  facility  location  models 
is  that  although  primarily  quantitatively  based,  qualitative 
factors  can  be  input  and  evaluated  in  some  formulations 
(10) .  Qualitative  factors  could  include  such  things  as  does 
a  site  have  favorable  tax  laws,  a  large  labor  pool,  or 
access  to  recreational  activities?  Researchers  speaking  on 
the  added  flexibility  of  location-allocation  models,  expound 
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on  the  fact  that  these  models  have  the  ability  to  represent 
wide  ranging  environments  in  mathematical  terms  (9). 

For  those  interested  in  a  brief  history  of  the  classic 
Euclidean  minisum  distance  facility  location  problem  or  more 
succinctly  known  as  the  Weber  problem,  Wesolowsky's  "The 
Weber  Problem:  History  and  Perspectives,"  is  a  good 
starting  point  (11) .  He  breaks  down  problem  development 
chronologically  and  credits  those  who  contributed  to  its 
present  form  and  understanding.  From  Fermat  (1601-1665), 
Torricelli  (1608-1647),  and  Cavalieri  (1598-1647)  to  the 
many  others  from  the  seventeenth  through  twentieth 
centuries,  Wesolowsky  gives  a  substantial  overview  of  the 
spatial  median  (Weber)  problem  history.  He  then  discusses 
the  generalization  of  the  Weber  problem  to  the  location- 
allocation  model,  where  points  to  be  located  are  facilities 
and  fixed  points  in  the  formulation  become  supply  and  demand 
points.  Transportation  costs  were  included  as  functions  of 
distance  and  mention  was  made  of  spherical  distances  as 
opposed  to  Euclidean  distances  (11) .  Continuing  with  the 
theme  of  location-allocation  models,  Ghosh  and  Harche  (9) 
also  review  their  progress  over  time. 

Ghosh  and  Harche  begin  from  the  introduction  of 
location-allocation  models  in  the  1960s  and  follow  their 
evolution  to  the  1990s.  They  cite  the  most  important 
characteristic  as  "the  ability  of  these  models  to  determine 
the  optimal  location  of  several  facilities  simultaneously." 
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In  many  distribution  systems  some  of  the  located  facilities 
are  used  as  transshipment  points  which  collect  goods  from 
dispersed  suppliers  and  then  ship  to  demand  points.  The 
same  objective  is  apparent  for  placing  these  transshipment 
centers  as  is  the  case  for  most  location  problems,  that 
being  the  minimization  of  cost  (9) .  Before  continuing  with 
transshipment  centers,  a  brief  overview  of  some  location 
analysis  that  has  been  applied  in  military  decision-making 
will  be  discussed. 

Despite  the  tremendous  amount  of  literature  on  facility 
location  models,  there  has  been  very  little  work  done  on  the 
placement  of  military  consolidation  points  within  the 
continental  United  States.  The  "Optimal  Airlift 
Distribution  Study  (1),"  discussed  in  Chapter  I,  was  the 
first  attempt  at  locating  CONUS  aerial  ports.  The  remaining 
studies  found  in  the  literature  search  that  concentrated  on 
military  location  applications  were  not  specifically 
directed  at  the  aerial  port  location  problem.  The  majority 
was  found  in  the  stack  of  theses  at  the  Air  Force  Institute 
of  Technology  (AFIT) .  Garcia  developed  and  applied  a 
coverage  type  location-allocation  problem  to  the  locating  of 
Air  Force  repair  facilities  and  the  associated  limited 
reparable  equipment  stocking  those  facilities  (12) .  Merrill 
tackled  a  single  facility  location  and  routing  problem  in 
order  to  site  and  minimize  the  en  route  distance  of  flight 
inspection  missions.  He  modified  and  applied  the  classic 
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"multiple  traveling  salesman"  model  and  solved  using  the 
Simplex  solution  method  (13)  . 

More  in  line  with  the  aerial  port  location  task  at 
hand,  a  thesis  was  done  on  the  location  and  routing  of  the 
Defense  Courier  Service  (DCS)  Aerial  Network.  The  DCS  is  an 
organization  whose  purpose  is  to  handle  and  transport 
sensitive  material  for  the  Department  of  Defense.  Again  the 
traveling  salesman  formulation  was  used  but  the  starting 
point  was  Laporte's  algorithm  input  as  an  integer  linear 
program  and  a  combined  heuristic  technique,  minimum  spanning 
forest/Clarke-Wright,  was  applied  to  obtain  a  solution  (14) . 

The  remaining  portion  of  this  chapter  will  briefly 
describe  two  formulations  of  the  multiple  transshipment 
center  location  problem.  These  models  could  be  applicable 
to  the  AMC  aerial  port  network  and  solved  within  a  given  set 
of  constraints. 

Model  1 ;  Automobile  Manufacturer 

In  1992  a  study  was  conducted  on  behalf  of  an  American 
automobile  manufacturer  to  locate  an  appropriate  number  of 
transshipment  centers.  These  centers  would  serve  as 
consolidation  points  for  small  Just-in-time  shipments  from 
hundreds  of  suppliers.  The  material  could  then  be 
transported  to  assembly  plants  in  a  more  cost-effective 
manner.  Bhaskaran' s  approach  to  this  particular  problem  was 
taken  as  a  continuous  space  model  as  opposed  to  a  network 
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model  (15) .  Network  models  have  a  specified  set  of 
alternate  locations  to  choose  the  best  from,  whereas 
continuous  models  can  locate  facilities  from  an  almost 
infinite  solution  space. 

Bhaskaran' s  objective  was  to  minimize  the  total  flow- 
weighted  transportation  distance  of  material  shipped  to  the 
plants.  The  formulation  of  this  large  problem  is  shown 
below : 


Objective  function: 


Minimize 

Z  -  wSkSiFkdkiSki  +  SjSiFjiDij 

(2.1) 

Subj  ect 

to : 

5ki  =  1  if  dki  <=  dki  for  all  1  9^  I,  and 

i  <=  1  for  all  1  such  that 
dki  —  dki, 

0  otherwise. 

(2.2) 

SiSki  =  1, 

(2.3) 

Fk  =  Ejfkj, 

(2.4) 

Fji  =  Skfkj  Ski 

(2.5) 

Where, 

oti,  pi 

=  location  of  center  I 

=  flow  from  supplier  k  to  customer  j 

dki 

=  distance  between  supplier  k  and  center 
(spherical  distance) 

I 

Dii 

=  distance  between  center  i  and  customer 

j 

w 

=  inbound  weight  factor 

Fk 

=  flow  from  supplier  k 
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=  flow  to  customer  j  through  center  I  (15) . 


Fji 


In  this  formulation,  w,  was  used  to  weight  inbound 
shipments  more  heavily  than  outbound  shipments.  Bhaskaran 
thought  this  appropriate  due  to  the  "[circuity]  of  inbound 
routes  and  the  loading  inefficiency  of  inbound  material 
(relative  to  outbound  material)  (15)." 

To  solve  this  problem,  Bhaskaran  used  a  multiple 
facility  heuristic  solution  procedure.  First,  for  a  given 
number  of  centers,  he  determined  the  best  locations.  That 
is  to  say  that  one-center,  two-center,  and  up  to  twenty- 
center  problems  were  solved.  Of  course  as  the  number  of 
centers  was  increased,  savings  in  total  ton-miles  were  seen, 
but  at  a  decreasing  rate.  Also,  centers  selected  from  early 
runs  remained  good  candidates  in  subsequent  runs  with  a 
greater  number  of  facilities  being  placed.  In  order  to 
choose  an  appropriate  final  number  of  transshipment  centers, 
he  introduced  a  minimum-size  requirement  and  reshuffled  the 
remaining  workload  as  the  smallest  centers  were  eliminated. 
Using  this  approximate  solution  method,  Bhaskaran  determined 
the  best  location  of  facilities  and  the  final  count  was  a 
total  of  eight  centers  (15) . 

His  approach,  as  previously  mentioned,  is  one  method  to 
solve  a  continuous  space  model.  In  the  case  of  aerial  port 
location  and  this  thesis,  a  finite  set  of  possible  APOEs  is 
given  and  the  best  locations  will  be  chosen  from  them. 
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Therefore  Bhaskaran's  transshipment  center  location 
formulation  and  solution  does  not  apply  to  this  discrete 
problem. 
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Model  2;  Multicoininodity  Distribution  System 


In  1974  a  paper  was  published  concerning  the  modeling 
of  a  very  complex  multiple  facility  location  problem.  This 
research  was  not  earth  shattering  but  was  a  different 
approach  as  compared  to  previously  applied  formulation 
techniques.  The  techniques  developed  by  Geoffrion  and 
Graves,  were  applied  and  gave  very  favorable  results  to 
large  scale  real  world  problems  (16) . 

Their  multicommodity  distribution  system  is  setup  to 
handle  a  large  number  of  commodity  types  produced  at  several 
plants.  The  goal  was  to  satisfy  known  customer  demand 
within  dispersed  zones  by  routing  the  shipment  of  various 
commodities  through  distribution  centers.  One  stipulation 
is  that  a  particular  customer  zone  be  assigned  to  one 
distribution  center.  By  consolidating  material  at  a  single 
facility,  economies  of  scale  can  be  realized  for  the  center 
to  customer  portion  of  the  shipment  (16) . 

The  problem  is  formulated  as  a  mixed  integer  linear 
program  and  is  illustrated  below: 

Objective  function: 

Minimi  ZOx  =>  O;  y,  z  =  0,1  Sljkl  CijkiXijj^i 

+  Ek  [fkZk  +  VkEii  Diiyki]  (2.6) 

Subject  to: 


Ekl 

^ijkl  Sij 

for 

all 

ij 

(2.7) 

I: 

Xijki  =  Diiyki 

for 

all 

ikl 

(2.8) 

Zk 

yki  =  1 

for 

all 

1 

(2.9) 
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for  all  k 


(2.10) 


YkZk  2il  DiiYki  <—  VkZk 

Linear  configuration  constraints  on  y  and/or  z. 

Where, 

i  =  coimnodity 

j  =  plant 

k  =  distribution  center  (DC)  sites 

1  =  customer  demand  zones 

Sij  =  supply  (production  capacity)  for  commodity 

i  at  plant  j 

Dii  =  demand  for  commodity  i  in  customer  zone  1 

Vk,  Vk  =  minimum,  maximum  throughput  for  a  DC  site 

fk  =  fixed  cost  of  DC  at  site  k 

Vk  =  variable  unit  cost  of  throughput  for  DC 

Cijki  =  average  unit  cost  of  producing  and  shipping 

commodity  from  plant  through  DC  to  zone  1 
Xijki  =  amount  of  commodity  shipped  from  plant 

through  DC  to  customer  zone  1 

Yki  =  1  if  DC  k  serves  1,  otherwise  0 

Zk  =  1  if  DC  is  acquired  at  k,  otherwise  0  (16) . 

The  significance  of  the  "ijkl"  subscript  variables 

according  to  the  authors  is  twofold.  First,  in  some 

applications  it  is  necessary  to  keep  track  of  where  the 

original  shipment  ended  up,  whereas  in  previous  models  the 

use  of  the  triple  subscript  lacked  this  flexibility.  Other 

models  used  separate  variables  for  plant  to  center  and 

center  to  customer  shipments  "linked  by  a  flow  conservation 

constraint."  The  second  reason  is  that  the  variables  make 

the  incorporation  of  direct  plant  to  customer  shipments  an 
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easy  matter  if  the  customer  does  not  also  receive  material 
from  a  distribution  center (16). 

The  overall  objective  was  to  meet  the  given  demands  of 
the  customer  at  the  least  total  distribution  cost  while 
satisfying  all  of  the  constraints.  A  discrete  set  of 
possible  locations  for  the  distribution  centers  was  given 
and  the  final  solution  is  a  subset  of  these,  with  particular 
sizes  of  facilities  solved  for  and  customer  zones  assigned 
to  them  exclusively. 

Summary 

A  basic  literature  review  was  conducted  and  reported 
within  this  chapter.  A  large  number  of  location  problems 
exist  in  the  literature  and  many  address  multiple 
transshipment  facility  location.  But  relatively  few  are 
applied  to  military  specific  examples.  This  is  not  a  major 
problem  because  existing  models  can  and  should  be  modified 
to  fit  any  number  of  real  life  situations. 

Of  the  models  investigated  in  this  literature  search, 
the  Geoffrion  and  Graves  formulation,  except  for  the 
multicommodity  count,  looks  like  the  best  formulation  for 
this  aerial  port  analysis.  The  next  chapter  will  address 
the  modifications  necessary  to  make  it  applicable  to  the 
aerial  port  location  study  and  the  data  required  as  input 
for  the  new  formulation. 
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Ill .  Methodology 


This  chapter  describes  the  methodology  used  in  the 
completion  of  this  aerial  port  location  analysis.  The 
aerial  portion  of  the  Defense  Transportation  System, 
consists  of  aerial  ports  of  embarkation,  aerial  ports  of 
debarkation,  and  final  consumption  locations  outside  the 
continental  US.  Feeding  cargo  to  the  APOEs,  are  the 
distribution  depots  and  other  supply  points  within  the 
CONUS.  In  this  study,  the  focus  is  on  that  portion  of  the 
system  composed  of  the  origination  cities,  APOE 
transshipment  bases,  and  APOD  arrival  points.  In  effect  the 
area  under  study  can  be  represented  as  a  distribution 
network  and  therefore  be  modeled  using  one  of  the  techniques 
explained  in  the  literature  review  of  Chapter  II. 

One  can  see  that  due  to  the  difficulty  and  exorbitant 
expense  of  establishing  or  moving  an  APOE,  there  would  only 
exist  a  select  few  locations  suitable  for  the  purpose.  Good 
candidates  for  basing  an  aerial  port  would  of  course  include 
the  facilities  already  established  by  the  DOD  along  the 
coastal  United  States.  Those  locations  are  Charleston  AFB, 
SC,  Dover  AFB,  DE,  McChord  AFB,  WA,  McGuire  AFB,  NJ,  Norfolk 
NAS,  VA,  and  Travis  AFB,  CA.  Three  additional  inland  sites 
were  chosen  to  include  in  the  formulation  as  alternatives  to 
the  current  structure.  Those  additional  sites  are  also 
established  Air  Force  bases  and  two  contain  Defense 
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Distribution  Depots  (cargo  origination  points) .  The  depot 
bases  are  Hill  AFB,  UT,  and  Tinker  AFB,  OK.  Two  of  the 
sites  also  contain  limited  aerial  port  facilities,  Tinker 
AFB,  OK  and  Wright-Patterson  AFB,  OH.  These  additional 
three  inland  sites  taken  with  the  six  established  APOEs, 
make  up  the  total  solution  set  considered  in  this  thesis. 

Cargo  shipment  data  was  taken  from  the  FY96 
Transportation  Reporting  and  Inquiry  System  database  and 
analyzed  to  incorporate  into  the  model.  The  summary  data 
was  inputted  to  represent  the  demand  at  APODs . 

Cost  and  distance  data  came  from  a  number  of  different 
sources.  They  will  be  discussed  separately  in  subsequent 
sections  of  this  chapter. 

Model  Formulation 

The  basic  structure  of  the  aerial  port  distribution 
system  is  closely  related  to  the  network  distribution  system 
modeled  by  Geoffrion  and  Graves.  Their  model  was  introduced 
in  Chapter  II.  This  multicommodity  distribution  system  is 
described  as  having  several  different  commodities  produced 
at  dispersed  locations.  The  commodities  are  shipped  via 
distribution  centers  to  satisfy  known  demands  within 
different  customer  zones.  Also  in  the  Geoffrion  and  Graves 
model,  the  stipulation  is  made  that  one  customer  location  is 
assigned  to  one  distribution  center.  This  would  allow  for 
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consolidation  of  material  and  therefore  favor  the 
realization  of  economies  of  scale  (16) . 

In  this  multicommodity  model,  possible  locations  of 
distribution  centers  are  given.  Operating  costs  for  each 
site  are  given  and  transportation  costs  are  assumed  to  vary 
linearly  with  distance  shipped.  The  overall  thrust  of  the 
problem  was  to  determine  the  number  and  location  of  the 
distribution  centers.  This  would  be  based  on  the  least-cost 
combination  of  establishing  facilities  and  shipping  cargo 
from  supply  point,  to  center,  to  customer  in  order  to  meet 
the  given  demand  (16) . 

Some  changes  to  the  model  per  se  need  to  be  made  in 
order  to  apply  to  the  military  APOE  location  problem.  First 
of  all,  the  multicommodity  aspect  can  be  likened  to  a 
particular  city  or  base  of  origin.  Cargo  from  one 
origination  point  is  supplied  to  various  APODs  based  on 
previously  shipped  quantities  of  material.  These  previous 
known  amounts  are  considered  a  type  of  "commodity"  supplied 
by  a  specific  origin  and  are  demands.  Second,  only  the  most 
significant  origination  and  destination  points  were  included 
in  the  analysis.  This  reduced  the  complexity  of  the  problem 
and  still  allowed  a  reasonable  representation  of  the 
complete  system.  The  mixed-integer  linear  programming 
formulation  of  the  aerial  port  distribution  system  can  be 
written  as  follows: 

Objective  function: 
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Minimize  Z 

—  Zjk  (Cjk*djk  + 

Vj)*Xjk  Zkl  Ckl*dkl*Xkl 

t  Zk  fk*Zk 

(3.1) 

Subject  to: 

Sk  Xjk 

<=  Sj 

for  all  j 

(3.2) 

Zj  Xjk 

=  Zl  Xkl 

for  all  k 

(3.3) 

Zk  Xkl 

=  Di 

for  all  1 

(3.4) 

Zk  ykl 

=  1 

for  all  1 

(3.5) 

Zj  Xjk 

<=  Zk*Mk 

for  all  k 

(3.6) 

Xjk 

<=  yjk*Mk 

for  all  j,  k 

(3.7) 

all  variables  are  nonnegative 

Where , 

j  =  cargo  origin/supply  point, 

k  =  aerial  port  of  embarkation  (transshipment) , 

1  =  aerial  port  of  debarkation  (demand  point) , 

Cjk  =  weighted  average  cost  per  ton-mile  of  shipping 
from  any  origin  j  to  any  APOE  k, 

Cki  =  weighted  average  cost  per  ton-mile  of  shipping 
from  any  APOE  k  to  any  APOD  1, 

djk  =  statute  mile  distance  from  origin  j  to  APOE  k, 

dki  =  nautical  mile  distance  from  APOE  k  to  APOD  1, 

Vj  =  APOE  throughput  cost  per  ton  of  cargo, 

Xjk  =  flow  in  tons  per  month  of  cargo  shipped  from 

origin  j  to  APOE  k, 

Xki  =  flow  in  tons  per  month  of  cargo  shipped  from 
APOE  k  to  APOD  1 

f^s.  =  monthly  operating  cost  for  APOE  k 

Zk  =  a  0  -  1  variable;  1  if  APOE  is  established  at  k, 

and  0  otherwise. 
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Yki  =  a  0  -  1  variable;  1  if  APOE  k  serves  APOD  1,  and 
0  otherwise, 

Sj  =  origin  supply  limitation, 

Di  =  demand  at  APOD  1,  and 

Mk  =  maximum  OCONUS  throughput  of  APOE  k  in  tons  per 
month . 

The  data  requirements  for  the  above  formulation  was  inputted 
into  a  Microsoft  Excel  (21)  spreadsheet  and  via  a  macro, 
output  in  a  format  for  the  CPLEX  (22)  linear  programming 
package  to  solve.  The  formulation  was  then  read  into  CPLEX. 
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Cargo  Data 


The  analysis  and  breakdown  of  cargo  demand  data  began 
with  the  FY96  TRAIS  database.  The  shipments  contained 
within  this  database  are  uniquely  identified  by  a  seventeen 
digit  transportation  control  number  (TCN) .  The  TCN  is  a  an 
alpha-numeric  code  used  by  the  DOD  in  accordance  with  the 
Military  Standards  and  Movement  Procedures  (MILSTAMP)  to 
identify  individual  shipments  within  the  DTS.  The  portion 
of  the  TRAIS  database  used  for  this  study  was  updated  on  27 
February  1997  and  stored  on  computer  disk  by  AMC  (4  and  5) . 

To  facilitate  extraction  of  useful  information  from 
TRAIS,  it  was  necessary  to  obtain  two  other  types  of 
identification  data.  They  were  Department  of  Defense 
Activity  Address  Codes  (DODAACs)  and  Air  Terminal 
Identification  Codes  (ATICs) .  The  DODAAC  "is  a  six  position 
alpha  numeric  code  that  identifies  a  unit,  activity,  or 
organization  that  has  the  authority  to  requisition  and/or 
receive  material  from  DOD."  (23).  The  DODAAC  is  further 
separated  into  types  of  address  codes  (TACs) :  TAC  1, 
identifies  a  unit's  mailing  address;  TAC  2,  is  a  freight 
address;  and  TAC  3,  is  a  billing  address,  which  in  many 
cases  is  at  a  location  hundreds  of  miles  away  from  the  cargo 
delivery  site.  This  separation  initially  caused  problems 
accurately  identifying  cities  of  origin.  The  problem  was 
resolved  with  the  help  of  the  organization  responsible  for 
maintaining  the  DODAAC  database,  the  Defense  Activity 
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Address  System  Center  (DAASC)  located  at  Wright-Patterson 
AFB,  OH  (21) .  ATIC  codes  were  provided  by  AMC  (5) .  An  ATIC 
is  the  three-letter  code  identifying  a  unique  airfield  in 
the  world. 

The  FY96  TRAIS  is  a  very  large  database  and  for  this 
analysis  began  with  almost  two  million  entries.  Microsoft 
Access  (22)  database  management  system  was  used  in 
conjunction  with  Microsoft  Excel  to  manage  the  data.  One 
quarter  of  cargo  data  was  extracted  from  the  original  table, 
as  this  would  be  used  to  represent  cargo  movement  for  the 
year.  Duplicate  TCNs  were  then  removed.  The  next  step  was 
to  limit  data  to  only  the  cargo  that  transited  the  six  major 
APOEs  of  Charleston,  Dover,  McChord,  McGuire,  Norfolk,  and 
Travis . 

Linking  the  resultant  table  containing  outbound  cargo 
data  with  both  the  DODAAC  information  and  ATIC  table, 
queries  were  run  to  find  the  biggest  shippers  and  ports  of 
debarkation  in  the  DTS  served  by  the  current  APOEs.  A 
consignor  list  by  DODAAC  was  generated  and  the  top  origins 
were  chosen  to  represent  approximately  85  percent  of  the 
weight  shipped  through  the  six  major  APOEs.  Consignors 
within  the  same  or  near  cities  were  consolidated,  as  this 
should  not  significantly  effect  the  final  locations  of  the 
APOEs  and  it  also  served  to  reduce  the  complexity  of  the 
task.  The  number  of  destinations  was  also  limited  to  the 
top  APODs  that  represented  approximately  85  percent  of  the 
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weight  delivered  outside  the  CONUS  via  the  six  major  APOEs. 

A  couple  of  possible  APOD  sites  were  eliminated  from  the 
analysis  based  on  blank  identification  locations  within  the 
TRAIS. 

Matching  consignor  DODAACs  with  APODs  and  consolidating 
cities  of  origin,  cargo  movement  information  between  these 
sites  was  generated  and  represented  approximately  75  percent 
of  the  weight  shipped  through  the  APOEs.  This  new  data  was 
adjusted  to  equal  100  percent  of  the  weight  transiting  the 
CONUS  APOEs  outbound.  The  new  figures  were  then  used  as 
demand  at  particular  APODs  from  specific  origins,  and  that 
demand  must  be  met  in  the  linear  program. 


25 


Distance  Data 


Distance  between  cargo  origin  points  and  APOEs  were 
taken  from  the  military  regulation.  Transportation  and 
Travel  Official  Table  of  Distances  (23) .  The  tabled  figures 
in  this  publication  are  based  on  driving  distance  and  are 
statute  miles. 

Distances  from  CONUS  APOEs  to  OCONUS  APODs  were  for  the 
most  part  obtained  from  a  Borland  Dbase  IV  database  used  at 
AMC  (24) .  The  database  was  converted  from  the  Dbase  IV  to 
Microsoft  Access  format  and  appropriate  distances  in 
nautical  miles  were  extracted.  In  some  cases  where  no  APOE- 
APOD  match  was  shown,  it  was  necessary  to  use  the  great 
circle  equation  to  calculate  distances.  The  formula  listed 
below  is  from  an  Air  Force  air  navigation  manual: 

D  =  60cos“^  [sin  (lati)  sin  (lat2) 

+  cos  (lati)  cos  (lata)  cos  (long2  -  longi)  ]  (3.8) 

The  lat  and  long  represent  the  latitude  and  longitude, 
respectively,  and  degrees  must  be  converted  to  radians  in 
order  to  use  this  formula.  The  distance,  D,  is  in  nautical 
miles  (nm)  (25) . 
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Cost  Data 


For  CONUS  truck  transportation  costs,  an  average  ton- 
mile  rate  was  obtained  using  the  Military  Traffic  Management 
Command's  (MTMC)  Draft  "Traffic  Management  Progress  Report" 
(TMPR)  for  FY95  (26) .  FY95  data  was  used  because  it  is  from 
the  most  current  report  available.  Also,  as  per  notes  on 
the  TMPR,  seven  months  of  data  was  unrecoverable  and  was  not 
used  in  MTMC's  calculation  of  costs. 

The  percentage  of  shipments  in  the  truckload  (TL) 

(10,000  pounds  and  over)  and  less-than-truckload  (LTL)  (less 
than  10,000  pounds)  categories  were  determined.  These 
percentages  were  used  along  with  the  average  ton-mile  rate 
per  weight-break  to  calculate  a  weighted  average  cost  per 
ton-mile . 

For  air  cargo  transportation  costs,  Defense  Business 
Operations  Fund  -  Transportation  (DBOF-T)  airlift  rates  were 
taken  from  the  "US  Government  Department  of  Defense  (DOD) 
Rate  Tariffs"  appendix  of  the  DBOF-T  rate  guide  (27) .  The 
rates  are  broken  down  into  five  weight-breaks  and  are  listed 
as  dollars  per  pound-mile.  These  rates  were  converted  into 
dollars  per  ton-mile  and  then  a  weighted  rate  per  ton-mile 
was  calculated.  The  breakdown  of  the  percentage  of 
shipments  in  each  weight-break  category  was  performed  by 
AMC's  Cargo  Movement  Branch  using  the  FY96  TRAIS  database 
(28)  . 


27 


Aerial  port  operating  costs  were  obtained  from  AMC's, 
Financial  Management  and  Budget  Directorate  (HQ  AMC/FMBT) . 
The  costs  include  FY96  operating  costs  on  file  for  the  six 
major  ports  and  "approximately  $1.2  million  related  to 
maintenance  and  repair  at  these  aerial  ports.  There  are  no 
military  costs  [personnel]  included  in  these  numbers" (29) . 
Because  this  thesis  was  to  examine  alternative  basing 
locations  for  the  CONUS  APOEs,  the  operating  costs  for  the 
three  inland  port  sites  are  not  Jcnown.  Therefore  random 
numbers  between  the  highest  and  lowest  operating  costs  of 
the  known  APOEs  were  generated  and  used  for  the  three 
additional  inland  port  sites. 

Summary 

The  demand  data  with  the  restricted  set  of  locations, 
along  with  the  distance  and  cost  information  just  discussed, 
will  be  applied  to  the  modified  distribution  center  location 
model.  This  model  most  closely  matches  the  current  setup  of 
the  aerial  port  system  and  should  provide  some  insight  into 
the  problem.  The  results  from  the  application  of  the 
aforementioned  methodology  will  be  discussed  in  Chapter  IV. 
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IV.  Results 


This  chapter  will  present  the  findings  as  discovered  by 
the  application  of  the  methods  discussed  in  Chapter  III. 

The  research  questions  of  Chapter  I  will  form  the  heart  of 
the  results.  Restated,  they  are: 

1.  Based  on  a  modified  distribution  system  location 
problem,  what  are  the  optimal  locations  and  how  many 
aerial  port  facilities  are  needed? 

2.  How  does  the  new  structure  compare  with  the  current 
APOE  structure  on  a  cost  basis? 

Cargo  Data  Analysis 

As  stated  in  the  methodology,  the  original  cargo  data 
provided  by  AMC  was  the  FY96  TRAIS  database  containing 
1,916,541  line  entries,  with  identifying  transportation 
control  numbers  and  other  accompanying  information. 

Although  the  TRAIS  database  was  last  updated  on  27  February 
1997,  entries  can  still  be  deleted  or  added  by  AMC  until  an 
official  close-out  date  is  established.  As  of  August  1997, 
one  had  not  yet  been  set  (4) .  The  data  is  therefore  not  as 
accurate  as  it  could  be,  but  that  should  have  little  if  any 
impact  on  the  results  of  this  analysis. 

The  breakdown  of  cargo  began  with  the  use  of  Microsoft 
Access  to  limit  the  large  volume  of  data  to  a  representative 
set  of  one  quarter  of  FY96.  That  quarter  was  arbitrarily 
chosen  to  run  from  1  April  to  30  June  1996.  This  reduced 
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the  number  of  entries  to  484,704,  or  as  expected  to 
approximately  25  percent.  Next,  in  order  to  avoid  duplicate 
information  and  provide  channel  summary  data,  line  entries 
with  duplicate  TCNs  were  removed.  This  again  resulted  in  a 
substantial  reduction  in  the  set  under  study  to  374,791,  or 
about  77  percent  of  that  quarter's  information. 

One  of  the  stipulations  of  this  aerial  port  analysis 
was  to  move  the  same  tonnage  of  cargo  that  transited  the  six 
major  APOEs  for  overseas  delivery.  Therefore  the  cargo  data 
set  was  further  restricted  to  those  shipments  that 
originated  within  the  CONUS  and  transited  the  six  major 
APOEs  destined  for  overseas  APODs .  Those  APOEs  again  are 
Charleston  AFB,  SC,  Dover  AFB,  DE,  Norfolk  NAS,  VA,  McChord 
AFB,  WA,  McGuire  AFB,  NJ,  and  Travis  AFB,  CA.  This  left  the 
total  number  of  line  entries  at  212,197. 

Using  the  relational  database  capabilities  of  Microsoft 
Access,  the  table  containing  the  212,197  entries  was  linked 
with  both  the  revised  DODAAC  table  and  the  ATIC  table. 
Queries  were  then  run  to  find  the  biggest  shippers  and  ports 
of  debarkation  by  weight.  The  first  query  resulted  in  a 
consignor  list  of  12,136  entries  separated  by  service,  APOE, 
APOD,  and  city.  The  next  consolidation  of  data  was 
accomplished  by  summing  all  of  the  cargo  originating  from 
the  same  consignor.  The  list  was  again  queried  to  show  the 
top  consignors  by  DODAAC,  and  the  top  93  are  shown  in 
Appendix  A.  The  cargo  generated  by  these  origins  represent 
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84.61  percent  of  the  weight  shipped  out  of  the  CONUS  through 
the  six  major  APOEs.  The  one-quarter  tonnage  shipped 
through  these  ports  was  23,614  for  a  monthly  average  of 
7,871  tons.  Those  93  origins  were  consolidated  based  on 
their  proximity  to  one  another  and  the  resulting  list  was 
reduced  to  the  53  shown  by  city  and  state  in  Table  1. 

The  aerial  ports  of  debarkation  are  shown  in  Appendix 
B,  with  the  top  22  shown  here  in  Table  2.  The  top  22  APODs 
represent  88.77  percent  of  the  total  cargo  weight  delivered 
outside  the  continental  US  that  transited  the  six  major 
APOEs . 

Cross  referencing  the  top  53  cities  of  origin  with  the 
top  22  APODs  using  Access,  resulted  in  the  extraction  of 
75.07  percent  of  the  cargo  which  originated  within  the  CONUS 
and  was  shipped  overseas  via  the  major  APOEs.  This  cargo 
information  by  weight  is  shown  in  Appendix  C.  In  order  to 
put  a  more  accurate  load  into  the  model  of  the  aerial  port 
system,  the  75  percent  tonnage  figures  were  increased  to 
equal  100  percent  of  the  cargo  originating  within  the  CONUS 
that  was  shipped  overseas  via  the  major  ports.  The  nodes 
were  not  changed  but  the  shipping  weights  were  modified. 

Distance  Data 

The  statute  miles  between  cities  of  origin  and  APOEs, 
taken  from  the  Transportation  and  Travel  Official  Table  of 
Distances,  are  shown  in  Table  3.  For  those  locations  not 
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included  in  the  regulation,  it  was  necessary  to  obtain  the 
distance  to  a  nearby  city  and  manually  adjust  the  mileage. 

Nautical  miles  between  APOEs  and  APODs,  as  taken  from 
an  AMC  table  of  distances,  are  shown  in  Table  4. 
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Table  1.  Top  Cities  of  Origin, 
(1  April  -  30  June  1996) 


CITY  OF  ORIGIN 


1  Anniston 


2  Ft  Rucker 


3  Huntsville 


Ft  Huachuca 


6  ILathrop 


Lemoore  NAS 


Monterrey 


5  lOakland 


10  San  Diego 


11  Travis  AFB 


12  Peterson  AFB 


27  Washington 


13  Dover  AFB 


14  EglinAFB 


15  Jacksonville 


16  Orlando 


17  Ft  Banning 


19  Ft  Stewart 


18  Palmetto 


20  Robins  AFB 


22  Rock  Island 


23  Scott  AFB 


24  Crane 


25  Ft  Campbell 


26  Ft  Knox 


28  Kessler  AFB 


29  Malmstrom 


30  Camp  Lejuene 


OffuttAFB 


34  I  McGuire  AFB 


35  Nellis  AFB 


33  New  York 


36  Columbus 


37  Tinker  AFB 


38  New  Cumberland 


39  IPhlladelphia 


40 


41  jCharleston 


42 


43  Memphis 


44  Corpus  Christ! 


45  Fort  Worth 


46  Ft  Hood 


47  San  Antonio 


48  Texarkana 


49  Hill  AFB 


51  Norfolk 


50  Richmond 


52 


53  Oak  Harbor 


STATE 


lAL 


AL 


AL 


AZ 


CA 


CA 


CA 


ICA 


CA 


CA 


CA 


CO 


pc  _ 


pE 


IfL _ 


FL 


FL 


GA 


GA 


GA 


GA 


IL 


|IL 


IL 


IN 


KY 


KY 


MS 


MT 


NC 


NC 


INE 


NJ 


NV 


NY 


OH 


OK 


PA 


PA 


PA 


ISC 


SC 


ITN 


TX 


TX 


TX 


TX 


TX 


UT 


VA 


VA 


WA 


jWA 


WEIGHT  (lbs. 


308,933 


182,473 


93,9461 


2,441,370 


84,201 


717,128 


141,391 


565,434 


794,713 


1,793,062 


102,938 


746,117 


2,420,973 


358,817 


114,313 


810,778 


135,608 


180,196 


257,477 


271,670 


158,156 


256,736 


260,601 


130,099 


81,383 


99,325 


71,730 


189,144 


286,984 


91,313 


1,987,133 


97,741 


183,029 


144,464 


377,168 


12,831,930 


196,056 


469,437 


1,244,857 


130,678 


112,748 


592,005 


143,238 


607,352 


781,398 


3,192,288 


1,052,708 


756,699 


75,883 


WEIGHT  (lbs.) 

39,959,720 


SHIPMENTS 


288 


236 


78 


241 


21,960 


282 


3,859 


117 


1,663 


4,952 


2,205 


242 


1,065 


4,718 


761 


268 


2,450 


267 


386 


818 


276 


328 


166 


212 


68 


339 


648 


131 


3,662 


222 


93 


6,499 


2,078 


54,313 


SHIPMENTS 

163,803 
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Table  2.  Top  Aerial  Ports  of  Debarkation 
(1  April  -  30  June  1996) 


APOD/ATIC 

CITY 

COUNTRY 

WEIGHT  (lbs.) 

SHIPMENTS 

1 

RMS 

RAMSTEIN  AB 

GERMANY 

36,528 

2 

OSN 

OSAN  AB 

SOUTH,  KOREA 

3,681,797 

8,986 

3 

HIK 

HONOLULU 

UNITED  STATES 

2,183,901 

10,408 

4 

OKO 

TOKYO 

JAPAN 

2,113,672 

15,206 

5 

KWI 

KUWAIT  CITY 

KUWAIT 

1,945,244 

1,269 

6 

HOW 

HOWARD  AB 

PANAMA 

1,567,047 

5,636 

7 

DHAHRAN 

SAUDI  ARABIA 

1,452,686 

11,102 

8 

KADENAAB 

JAPAN 

1,261,717 

10,369 

SIZ 

SIGONELLAAB 

ITALY 

1,180,224 

15,192 

MILDENHALLAB 

UNITED  KINGDOM 

1,060,807 

9,645 

11 

KEF 

KEFLAVIK 

ICELAND 

1,033,173 

3,681 

m 

BAH 

BAHRAIN 

BAHRAIN 

992,246 

7,851 

13 

NBW 

GUANTANAMO  BAY 

CUBA 

737,502 

1,494 

14 

NRR 

ROOSEVELT  ROADS 

PUERTO  RICO 

m 

ANDERSON  AFB 

GUAM 

629,812 

6,112 

THU 

THULE  AB 

GREENLAND 

579,544 

1,194 

RTA 

ROTA  (NAS) 

SPAIN 

5,657 

EDF 

ANCHORAGE 

UNITED  STATES 

472,560 

1,992 

19 

RUH 

_ 

RIYADH 

SAUDI  ARABIA 

432,963 

3,728 

20 

KWA 

KWAJALEIN 

US  TERRITORY 

398,640 

474 

m 

PLA 

PALMEROLA 

HONDURAS 

389,241 

415 

ASP 

ALICE  SPRINGS 

AUSTRALIA 

378,266 

131 

WEIGHT 

SHIPMENTS 

TOTAL 

HHHH 

Those  origin-destination  pairs  not  listed  in  AMC' s  mileage 
table  were  calculated  from  the  great  circle  equation  (3.8) 
and  are  shown  with  an  asterisk. 


34 


Table  3.  Origin  to  APOE  Distance 


CHS 

DOV 

HIF 

NGU 

TCM 

WRI 

SUU 

TIK 

FFO 

1 

Anniston  AL 

375 

813 

1850 

662 

2605 

879 

2392 

750 

547 

2 

Palmetto  GA 

304 

738 

1923 

587 

2678 

804 

2488 

846 

558 

3 

Benning  Ft  GA 

379 

840 

1957 

689 

2712 

906 

2492 

857 

657 

4 

Bragg  Ft  NC 

206 

426 

2162 

242 

2854 

496 

2814 

1172 

552 

5 

Campbell  Ft  KY 

587 

817 

1611 

733 

2366 

880 

2303 

661 

346 

6 

Charleston  SC 

10 

603 

2203 

416 

2955 

675 

2768 

1126 

697 

7 

Chicago  IL 

907 

779 

1406 

877 

2042 

778 

2111 

795 

286 

8 

Columbus  OH 

658 

505 

1699 

583 

2356 

504 

2404 

937 

67 

9 

Corpus  Christ!  TX 

1292 

1725 

1507 

1590 

2332 

1788 

1876 

594 

1334 

10 

Crane  IN 

644 

659 

1476 

697 

2221 

581 

2183 

687 

208 

11 

Dover  AFB  DE 

609 

10 

2162 

198 

2819 

111 

2867 

1417 

569 

12 

Eglin  AFB  FL 

506 

1036 

1980 

883 

2786 

1102 

2444 

877 

767 

13 

Fort  Worth  TX 

1099 

1455 

1253 

1371 

2078 

1518 

1717 

210 

1024 

14 

Harrisburg 

621 

146 

2028 

294 

2685 

134 

2733 

1305 

236 

15 

Hill  AFB  UT 

2195 

2162 

10 

2227 

825 

2161 

723 

1152 

1615 

16 

Hood  Ft  TX 

1169 

1572 

1306 

1448 

2121 

1635 

1748 

354 

1141 

17 

Huachuca  Ft  AZ 

1983 

2334 

874 

2250 

1593 

2358 

939 

930 

1795 

18 

Huntsville  AL 

470 

804 

1767 

710 

2522 

8  67 

2309 

667 

447 

19 

Jacksonville  FL 

262 

826 

2221 

627 

2976 

885 

2734 

1124 

867 

20 

Kessler  AFB  MS 

655 

1116 

1852 

965 

2671 

1182 

2316 

749 

835 

21 

Knox  Ft  KY 

641 

733 

1638 

669 

2369 

736 

2337 

791 

192 

22 

Lejuene  Camp  NC 

234 

389 

2288 

202 

2976 

479 

2936 

130 

674 

23 

Lemoore  NAS  CA 

2758 

2877 

833 

2793 

919 

2887 

209 

1452 

2324 

24 

Malmstrom  MT 

2279 

2151 

553 

2249 

668 

2150 

1120 

1402 

1658 

25 

McChord  AFB  WA 

2947 

2819 

825 

2917 

10 

2818 

724 

1968 

2326 

26 

McClellan  AFB  CA 

2759 

2815 

671 

2880 

721 

2814 

54 

1641 

2306 

27 

McGuire  AFB  NJ 

680 

111 

2161 

288 

2818 

10 

2866 

1438 

568 

28 

Memphis  TN 

665 

1  964 

1553 

!  880 

2308 

1027 

2095 

453 

533 

29 

Nellis  AFB  NV 

2256 

2510 

436 

2471 

1127 

2509 

569 

1138 

1964 

30 

New  York  NY 

735 

166 

2180 

346 

2828 

69 

2885 

1473 

603 

31 

Norfolk  VA 

422 

198 

2227 

10 

2917 

288 

2932 

1333 

617 

32 

Oakland  CA 

2752 

2911 

767 

2965 

768 

2910 

47 

1632 

2402 

33 

Offutt  AFB  NE 

1288 

1235 

930 

1300 

1680 

1234 

1635 

447 

726 

34 

Ord  Ft  CA 

2719 

3014 

880 

2926 

874 

3013 

153 

1593 

2465 

35 

Orlando  FL 

‘  392 

956 

2334 

769 

3093 

1027 

1  2828 

1234 

988 

36 

Peterson  AFB  CO 

1676 

1762 

580 

1729 

1396 

1761 

1269 

606 

1198 

37 

Philadelphia  PA 

650 

81 

2133 

271 

2790 

34 

2838 

1410 

540 

38 

Richmond  VA 

405 

205 

2139 

88 

2829 

275 

2844 

1256 

529 

39 

Robins  AFB  GA 

292 

791 

2006 

624 

2761 

857 

2571 

929 

641 

40 

Rock  Island  IL 

1038 

928 

1240 

1011 

1940 

927 

1945 

714 

435 

41 

Rucker  Ft  AL 

439 

960 

1968 

801 

2722 

1026 

2477 

875 

702 

42 

San  Antonio  TX 

1266 

1693 

1356 

1576 

2181 

1756 

1725 

481 

1262 

43 

San  Diego  CA 

2410 

2761 

775 

2677 

1217 

2785 

509 

1357 

2222 

44 

Scott  AFB  IL 

827 

923 

1376 

880 

2131 

922 

2075 

529 

370 

45 

Shaw  AFB  SC 

107 

558 

2137 

371 

2884 

629 

2729 

1087 

62  6 

46 

Stewart  Ft  GA 

154 

718 

2148 

523 

2902 

789 

2712 

1078 

739 

47 

Lathrop  CA 

2708 

2868 

724 

2923 

762 

2867 

52 

1590 

2359 

48 

Red  River  TX 

925 

1259 

1395 

1169 

2222 

1322 

1861 

288 

828 

49 

Tinker  AFB  OK 

1118 

1417 

1152 

1333 

1968 

1438 

1642 

10 

875 

50 

Tobyhanna  PA 

744 

184 

2101 

371 

2748 

126 

2805 

1414 

546 

51 

Travis  AFB  CA 

2760 

2867 

723 

2932 

724 

2866 

10 

1642 

2358 

52 

Washington  D.C. 

512 

103 

2072 

191 

2728 

169 

2776 

1314 

478 

53 

Whidbey  Island  WA 

2985 

2857 

862 

2955 

108 

2856 

824 

2007 

2364 

Source:  AFR  177-135  Transportation  and  Travel  Official 
Distances  are  in  statute  miles.  CHS-Charleston  AFB  SC, 
NGU-Norfolk  NAS  VA,  TCM-McChord  AFB  WA,  WRI -McGuire  AFB 
AFB  OK,  FFO-  Wright-Patterson  AFB  OH. 


Table  of  Distances. 

DOV-Dover  AFB  DE,  HIF-Hill  AFB  UT, 
NJ,  SUU-Travis  AFB  CA,  TIK-Tinker 
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Table  4 .  APOE  to  APOD  Distance 


APOD  APOE 


CHS 

DOV 

HILL 

NGU 

TCM 

WRI 

SUU 

TIK 

FFO 

1 

ASP 

8947 

*9063 

*7413 

9070 

6977 

*9078 

7585 

*8083 

*8672 

2 

BAH 

6340 

6391 

6659 

6291 

6396 

5948 

8593 

6695 

*6107 

3 

DHA 

6268 

5955 

6653 

6092 

6393 

5890 

7395 

6893 

6097 

4 

DNA 

7367 

6797 

5581 

6784 

5003 

7495 

5309 

6318 

6466 

5 

EDF 

3369 

2958 

1827 

3040 

1266 

3363 

1799 

2506 

2679 

6 

HIK 

4173 

4274 

2606 

4265 

2310 

4301 

2115 

3262 

3874 

7 

HOW 

1641 

1868 

2582 

1892 

3162 

1880 

3492 

2310 

1864 

8 

KEF 

2806 

2371 

3222 

2503 

4402 

2305 

4491 

3421 

2562 

9 

KWA 

*6236 

6400 

4626 

6298 

4774 

*6318 

4241 

5388 

*5916 

10 

KWI 

6059 

6167 

6451 

5757 

6419 

5573 

7798 

7195 

5885 

11 

MHZ 

3552 

3134 

4203 

3856 

5166 

3069 

5244 

4032 

3402 

12 

NBW 

820 

1158 

*2262 

1021 

*2823 

1826 

3383 

2201 

1196 

13 

NRR 

1171 

1389 

2738 

1252 

3266 

1672 

3229 

3297 

*1605 

14 

OKO 

7070 

5978 

4763 

6115 

4184 

6676 

4490 

5499 

5688 

15 

OSN 

7149 

6579 

5116 

6191 

4785 

7277 

5091 

6100 

5896 

16 

PLA 

1195 

1614 

2074 

1479 

2623 

1680 

2314 

1364 

1533 

17 

RMS 

3938 

3449 

4519 

3564 

5494 

3384 

4914 

4355 

3722 

18 

RTA 

3556 

3500 

4617 

3288 

4732 

3148 

5614 

4750 

3557 

19 

RUH 

6581 

5800 

6693 

5914 

6466 

6099 

8820 

7192 

6066 

20 

SIZ 

4577 

4446 

*5250 

4309 

6253 

4194 

5720 

5218 

4388 

21 

THU 

2702 

2276 

2400 

2389 

4298 

2201 

4387 

3317 

2244 

22 

UAM 

7468 

8992 

5451 

6993 

5544 

7596 

5410 

6557 

6605 

Source:  AMC  Mileage  Table  (Borland  Dbase  IV  File  converted  to  Microsoft  Access  95). 
Distances  are  in  nautical  miles.  Note:  *From  the  great  circle  equation. 

For  APOE  abbreviations  see  Table  3. 

ASP-Alice  Springs  Australia,  DAH-Dhahran  Saudi  Arabia, 

DNA-Kadena  AB  Japan,  EDF-Elmendorf  AFB  AK,  HIK-Hickam  AFB  HI 
HOW-Howard  AFB  Panama,  KEF-Keflavik  Iceland,  KWA-Kwajalein  Marshall  Island, 
KWI-Kuwait  City  Kuwait,  MHZ-Mildenhall  AFB  England,  NBW-Guantanamo  Bay  Cuba,  NRR- 
Roosevelt  Roads  NAS  Puerto  Rico,  OKO-Yokota  AFB  Japan,  OSN-Osan  AB  Korea,  PLA-Soto 
Cano  Honduras,  RMS-Ramstein  AB  Germany,  RTA~Rota  NAS  Spain,  RUH-Riyadh  Saudi 
Arabia,  SIZ-Sigonella  Italy,  THU-Thule  AB  Greenland,  UAM~Andersen  AFB  Guam. 


Cost  Data 


The  inland  truck  freight  transportation  costs  obtained 
from  MTMC's  "Traffic  Management  Progress  Report"  are  shown 
in  Table  5.  The  resultant  weighted  average  cost  was  0.2799 
dollars  per  ton-mile.  This  was  used  as  a  linear 
transportation  cost  function  for  shipping  from  city  of 
origin  to  aerial  port  of  embarkation. 
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Table  5.  Ton-mile  Costs  for  Inland  Freight  Traffic. 


Weighted 


Shipments  Average  Average 


TL  10,000  lbs  and  over 
LTL  Less  than  10,000  lbs 


Total 


(000)  Percent  Cost/Ton-mile  Cost/Ton-mile 


86.85 

19.03 

$0.0907 

369.54 

80.97 

$0.3244 

456.39 

100.00 

$0.2799 

Source:  MTMC's  Traffic  Manangement  Progress  Report  (Draft)  for  FY95. 

Data  from  March-September  excluded  due  to  reporting  deficiencies. 


Air  transportation  costs  are  also  dependent  on  the 
weight  break  category  that  a  shipment  falls  into.  The  AMC 
weight  breaks  used  in  this  problem  were  for  FY96  and  are 
shown  in  Table  6.  Also  shown  are  the  percentage  of 
shipments  within  each  weight  break.  One  can  see  that  the 
majority  of  the  shipments,  approximately  81  percent,  are  in 
the  smallest  category  of  1  to  439  pounds.  This  information 


Table  6.  AMC  Channel  Weight  Breaks  and  Costs 

Weighted 


Shipment  Size 
(lbs) 

Percent 

Rate 

Cost/Ton-mlle 

Average 

Cost/Ton-mlle 

1-439 

81.46 

440-1099 

10.75 

$0.9139566000 

1100-2199 

4.54 

$0.8139292000 

2200-3599 

1.48 

$0.7103598000 

3600  and  over 

1.77 

$0.6254112000 

TOTAL 

100 

$0.9846602192 

Source:  US  Government  -  DBOF-T  Airlift  Rate  Guide  FY96, 


is  again  shown  in  Figure  1  and  can  be  used  to  calculate  a 
weighted  average  of  air  transportation  cost.  The  resulting 


37 


value  of  0.9845  dollars  per  ton-mile  was  used  in  this 
application . 

Figure  1.  Shipments  within  Channel  Weight  Breaks 


Weight  Group  (lbs) 
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Source:  HQ/AMC/DONCM;  FY96  World-Wide  Channel  Shipment  Profile  Study 


Aerial  port  operating  costs  obtained  from  HQ  AMC/FMBT, 
are  shown  in  Table  7.  Values  drawn  randomly  from  the 
uniform  distribution  formed  by  the  highest  and  lowest 
port  operating  costs  at  Travis  and  McChord  were  used  as 
operating  costs  for  the  three  inland  bases  of  Hill, 
Tinker,  and  Wright-Patterson . 
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Table  7.  Aerial  Port  Operating  Costs  (FY96) . 


APOE 

Charleston  AFB  SC 
Dover  AFB  DE 
McChord  AFB  WA 
McGuire  AFB  NJ 
Norfolk  NAS  VA 
Travis  AFB  CA 
Hill  AFB  UT 
Tinker  AFB  OK 
Wright-Patterson  AFB  OH 


Operating  Cost  (000) 
$6,510.70 
$7,833.20 
$4,785.70 
$4,829.10 
$7,679.50 
$9,067.90 
$7,706.60 
$4,960.60 
$6, 613.70 


Source:  HQAMC/FMBT.  Note:  Hill,  Tinker,  and  Wright-Patterson 

costs  randomly  generated  from  uniform  distribution. 


Port  Capacities 

CONUS  aerial  port  throughput  capacities  shown  in  Table 
8,  are  all  based  on  the  current  manpower  authorized  (except 
Hill,  Tinker,  and  Wright-Patterson)  at  those  locations. 

That  is  to  say  that  if  manning  was  increased  during 
peacetime,  the  throughput  capability  would  also  increase. 
Also,  high  manning  levels  is  not  a  cure  all.  Other  factors 
such  as  material  handling  equipment,  ramp  space,  storage 
facilities,  and  fuel,  can  also  and  do  limit  the  capacity  of 
aerial  ports.  The  results  of  this  study  must  be  carefully 
examined  and  weighed  against  other  pertinent  variables 
before  any  final  decisions  are  made  on  the  future  of  the 
aerial  ports  of  embarkation  (30  and  31)  . 

The  total  throughput  capacity  figures  obtained  from  AMC 
are  shown  in  Table  8.  They  represent  the  amount  of  cargo  in 
tons  per  month  that  transit  these  aerial  ports  both  for 
CONUS  delivery  and  OCONUS  delivery.  Because  this  thesis  is 
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concerned  with  the  placement  of  CONUS  APOEs  not  APODs,  the 
total  capacity  of  these  ports  had  to  be  adjusted.  That 
adjustment,  results  shown  in  Table  8,  is  a  first  attempt  to 
capture  that  portion  of  a  port's  capacity  that  is  consumed 
by  the  outward  movement  of  cargo. 

As  with  port  costs,  capacities  of  the  three  inland 
bases  were  not  available  but  had  to  be  derived.  Drawing 
randomly  from  a  uniform  distribution  gave  the  results  in 
Table  8  for  Hill,  Tinker,  and  Wright-Patterson . 


Table  8.  Aerial  Port  Throughput  Capacities 


Conus  Aerial  Port 

%  as  APOE 

Total 

Throughput 

Capacitv 

(tons/month) 

OCONUS 

Throughput 

Caoacitv 

(tons/month) 

Charleston  AFB  SC 

70.60% 

5,500 

3,883 

Dover  AFB  DE 

67.86% 

9,500 

6,447 

Norfolk  NAS  VA 

66.78% 

6,000 

4,007 

McChord  AFB  WA 

67.81% 

1,500 

1,017 

McGuire  AFB  NJ 

81.08% 

2,000 

1,622 

Travis  AFB  CA 

58.26% 

8,500 

4,952 

Hill  AFB  UT 

4583 

Tinker  AFB  OK 

4363 

Wright-Patterson  AFB  OH 

3264 

Source:  HQ  AMC/DOZX  and  AMC  Air  Terminal  Norfolk  NAS. 


Hill,  Tinker,  and  Wright-Patterson  random  from  uniform 
distribution. 


40 


Old  vs.  New  Port  Structure 


The  detailed  data  discussed  throughout  this  chapter  was 
inputted  into  the  revised  formulation  of  a  distribution 
problem.  The  macro  written  by  Maj  Ray  Hill  which  creates  a 
CPLEX  readable  file,  is  shown  in  Appendix  D  (32) .  The 
results  of  the  subroutines,  those  being  the  objective 
function  and  constraints  of  the  linear  program  fed  into  the 
CPLEX  Linear  Optimizer,  can  be  provided  by  the  author  on 
request.  This  is  also  true  for  the  large  output  file  from 
CPLEX. 

The  first  run  of  the  formulation  indicates  that  three 
CONUS  aerial  port  facilities  should  remain  open,  two  on  the 
East  coast  and  one  on  the  West.  They  are  the  ports  at  Dover 
AFB,  DE,  McGuire  AFB,  NJ,  and  Travis  AFB,  CA.  The  total 
cost  calculated  for  that  system  was  33.41  million  dollars 
with  the  APOEs  sharing  the  monthly  workload  as  follows: 

Dover  -  4,212  tons,  McGuire  -  1,622  tons,  and  Travis  -  2,036 
tons.  The  total  cost  includes  CONUS  freight,  air  freight, 
and  port  operating  costs.  Naturally,  Travis  was  entirely 
focused  on  destinations  west  and  Dover  and  McGuire  to  the 
east.  Upon  examination  of  the  results,  one  notices  that 
McGuire  has  hit  its  capacity  limitation.  This  would  lead 
one  to  conclude  that  the  LP  opened  McGuire  first  because 
McGuire  had  the  least  expensive  operating  cost.  McGuire 
being  quickly  overwhelmed  with  cargo,  lead  to  the  opening  of 
Dover.  Intuitively  one  would  say  that  Dover  should  be  the 
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only  east  coast  port,  as  it  dwarfs  the  capacity  of  McGuire, 
could  handle  all  eastbound  cargo,  and  would  surely  limit  the 
overall  system  cost  if  McGuire  were  closed. 

In  order  to  test  the  above  theory,  the  LP  was  again 
solved,  only  this  time  Dover  was  forced  to  open  first.  The 
answer  returned  was  not  what  was  expected.  Dover,  Travis, 
and  McGuire  all  remained  open,  with  no  decrease  in  overall 
cost.  To  take  it  one  step  further  and  check  the  validity  of 
the  LP,  Dover  was  again  forced  open  but  this  time  McGuire 
was  forced  closed.  In  that  case  only  Dover  and  Travis  had 
the  active  aerial  ports,  but  overall  cost  increased. 

Although  a  relatively  minor  adjustment  of  about  49,000 
dollars,  the  increase  was  not  expected.  With  5,834  tons  and 
2,036  tons  transiting  Dover  and  Travis  respectively,  the 
system  cost  increased  to  33.46  million  dollars. 

The  decision  to  restructure  the  APOE  system  will  and 
should  be  based  on  many  factors.  For  this  thesis,  the  key 
factor  of  comparison  between  the  current  and  new  structure 
is  cost.  Those  are  port  operating  and  transportation  costs, 
both  CONUS  truck  freight,  and  overseas  air  freight.  The 
current  APOE  structure  was  compared  with  the  model  output 
based  on  the  cargo  data  set  used  for  this  thesis.  That 
would  include  7,871  tons  per  month,  port  operating  costs 
provided  by  AMC/FMBT,  0.2799  dollars  per  ton-mile  for  CONUS 
freight,  and  0.9845  dollars  per  ton-mile  for  air  freight. 
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The  current  DTS  structure  contains  the  six  coastal 
aerial  port  facilities  of  Charleston,  Dover,  McChord, 
McGuire,  Norfolk,  and  Travis.  Using  the  above  comparison 
factors,  a  rough  cost  estimate  was  obtained  for  today's 
system.  Those  costs  are  listed  in  Table  9.  Also  shown 
is  the  total  cost  of  the  three-port  structure.  The 
difference  in  monthly  cost  of  the  two  systems  is  960 
thousand  dollars  per  month  for  an  annual  figure  of  11.5 
million  dollars 

Table  9.  Old  vs.  New  APOE  Structure  -  Cost  Comparison 

System 

Six-Port  Three-Port 

Cost /Month 

Total  Freight  $30,979,421  $31,600,480 

Operating  $3,392,175  $1,810,850 

TOTAL  $34,371,596  $33,411,330 
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Summary 


Applying  the  methodology  from  Chapter  III  led  to 
answers  to  the  basic  questions  of  this  research  effort. 

Those  were: 

1.  What  are  the  optimal  locations  and  how  many  aerial 

port  facilities  are  needed? 

2.  How  does  the  new  structure  compare  with  the  current 

APOE  structure  on  the  basis  of  cost? 

By  analyzing  the  flow  of  cargo  from  CONUS  origination 
sites,  to  aerial  ports  of  debarkation,  the  data  points  input 
into  the  model  were  reduced  to  a  easily  manageable  number. 
The  final  cut  represented  over  75  percent  of  the  total  flow 
of  cargo  through  the  APOEs.  Distances,  costs,  and 
capacities  were  obtained  from  official  sources  and  when 
necessary  derived.  This  string  of  data  was  then  input  into 
the  modified  multicommodity  distribution  formulation  of 
Geoffrion  and  Graves,  leading  to  the  recommendation  to  only 
keep  the  three  APOEs  of  Travis,  Dover,  and  McGuire  open. 

The  model  portrays  an  approximate  one  million  dollar 
saving  if  the  new  structure  replaces  the  current  one.  Some 
of  the  other  factors  influencing  the  location  of  aerial 
ports  will  be  briefly  mentioned  in  conclusion. 
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V.  Conclusion 


With  the  continued  reduction  in  force  size  and  base 
structure,  the  US  Armed  Forces  retreat  more  and  more  to  the 
continental  United  States.  Exacerbating  the  situation  is 
the  relentless  budget  cutting  of  Congress.  As  a  result,  the 
requirement  for  the  current  high  level  of  material  support 
provided  to  overseas  locations  is  waning.  The  force 
restructuring  to  include  the  consolidation  of  distribution 
depots  under  one  "roof"  has  forced  another  look  at  the 
aerial  port  system.  More  specifically  where  and  how  many 
CONUS  aerial  port  facilities  should  the  future  Defense 
Transportation  System  have?  Once  that  was  determined,  how 
does  the  cost  of  the  new  system  compare  to  that  of  the  old? 
That  was  the  basic  thrust  of  this  thesis  and  answers  to  the 
above  questions  were  obtained  by  modeling  the  APOE  system  as 
a  transshipment  problem. 

This  effort  began  with  a  review  of  the  literature 
available  dealing  with  facility  location  problems.  Through 
extensive  research,  the  focus  was  narrowed  to  one  of  two 
possible  solution  methods,  or  models.  One,  by  Bhaskaran, 
was  applied  to  an  automobile  manufacture's  location  of 
transshipment  facilities  between  suppliers  and  final 
assembly  plants.  The  multicommodity  distribution  system 
formulation  of  Geoffrion  and  Graves  was  the  other  possible 
model.  It  was  determined  that  the  multicommodity  model 
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could  be  modified  to  best  represent  the  current  aerial  port 
system. 

The  result  of  this  modified  formulation  was  that  three 
CONUS  aerial  port  facilities  should  remain  open,  two  in  the 
East  and  one  in  the  West.  This  port  structure  allows 
increased  consolidation  of  cargo  for  air  shipment  and 
significantly  reduces  the  total  cost  of  shipment  as  defined 
in  this  problem.  The  savings  stem  from  the  operating  cost 
reduction  with  the  closure  of  the  three  APOEs  of  Charleston, 
McChord,  and  Norfolk. 

There  are  a  couple  of  other  benefits  in  this  three  port 
solution  than  just  cost.  With  the  opening  of  three  versus 
two  ports  (2  could  handle  the  capacity)  there  is  a  greater 
surge  capability  built  into  the  system.  Of  course  it  would 
not  be  as  high  as  a  six-port  system  would  but  the  capability 
to  handle  any  number  of  emergencies  or  contingencies  is 
still  there.  With  the  overall  cost  of  the  system  declining, 
AMC  could  move  for  a  reduction  in  the  price  of  air  shipment 
for  its  DOD  customers. 

It  must  be  realized  that  numerous  assumptions  and 
simplistic  representations  were  made  in  order  to  reduce  the 
complexity  of  the  aerial  port  location  study.  The  short 
list  includes  the  following.  The  number  of  nodes,  origins 
and  APODs,  were  limited  to  only  represent  75  percent  of  the 
tonnage  transiting  the  six  major  APOEs.  The  demands  in  the 
model  were  deterministic  and  were  exactly  met,  whereas  in 


46 


reality  there  is  always  some  variability.  The  only  port 
costs  considered  were  operating  costs.  They  did  not  include 
personnel  and  numerous  other  costs.  Capacities  used  were 
derived  based  on  the  percentage  of  outbound  cargo  a  port 
handles.  This  short  list  plus  other  limiting  factors 
omitted,  such  as  storage  capacity,  ramp  space,  fuel,  cargo 
handling  equipment,  and  other  support  activities,  must  be 
examined  carefully  before  any  closure  actions  are  taken. 
Major  areas  not  addressed  in  this  research  also  include 
aircraft  inventory  and  basing  strategies  and  route  timing 
and  location.  These  may  have  an  impact  on  CONUS  port 
locations,  and  would  definitely  need  to  be  addressed  in  a 
larger  more  comprehensive  study  of  the  whole  distribution 
system. 
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Appendix  A:  Top  Consignor  DODAACs . 
(1  April  -  30  June  1997) 


STATE  I  WEIGHT  (lbs.)  |  SHIPMENTS 


2  FB4497 


3  FB4427 


5  ISW3225 


SW3100 


8  SB3106 


10  SW0400 


11  S1002A 


12  SW3200 


13  FB4479 


14  SW3400 


15  W73BFY 


17  GN0003 


18  SW3224 


19  SW3213 


23  SW3218 


24  SW3114 


26  SW3500 


27  N00188 


28  SW3212 


29  SW3227 


39  FB4800 


M3 1000 


41  GA0001 


42  N65580 


New  Cumberland 


Dover  AFB 


Travis  AFB 


McGuire  AFB 


Lathrop 


Charleston 


Mechanicsburg 


Norfolk 


New  Cumberland 


Richmond 


Orlando 


Tracy 


McChord  AFB 


Ogden 


Arlington 


Norfoik 


Burlington 


Lathrop 


Kelly  AFB 


Norfolk 


McClellan  AFB 


Tinker  AFB 


San  Diego 


Tobyhanna 


Newport  News 


Memphis 


Norfolk 


McClellan  AFB 


Texarkana 


San  Diego 


Ft  Bragg 


Robins  AFB 


Mineral  Wells 


Crane 


Hill  AFB 


Langley  AFB 


[Camp  Lejuene 


Palmetto 


Portsmouth 


10,898,2221 


2,420,973 


1,793,0621 


1,407,893 


1,224,345 


1,058,155 


1,029,246 


982,940 


904,462 


611,161 


577,747 


564,546 


531,131 


487,944 


465,528 


462,378 


457,757 


456,967 


401,924 


377,168 


361,681 


355,286 


354,862 


334,053 


323,818 


315,204 


308,530 


295,772 


271,670 


268,536 


260,601 


35 

SCGA08 

Union  City 

CA 

258,410 

36 

W52H1C 

Rock  Island 

IL 

256,736 

37 

HXYAAA 

Dallas 

TX 

247,963 

216,852 


204,656 


189,144 


183,114 


3,372 


4,718 


2,205 


1,079 


2,007 


910 


24,264 


1,221 


26,677 


246 


3,919 


713 


2,039 


2,516 


7,255 


2,552 


7,143 


150 


2,078 


3,484 


1,061 


722 


3,814 


789 


3,709 


929 


1,306 
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43  IW31R4Z 


44  W33GGZ 


45  W31G1Z 


46  FB2823 


47  S1403A 


48  S4404A 


49  SW3120 


50  SW0700 


52  W45QRE 


53  W62R65 


55  W33APT 


56  FB4803 


57  W34GNA 


58  N57012 


59  SP0200 


60  FB4819 


61  W23A9F 


65  N60478 


66  SW3222 


67  N00228 


69  W68P4L 


70  SP0400 


72  N65236 


73  CL0KX3 


74  FB3010 


75  S2101A 


76  N44399 


77  FB4852 


78  W31P38 


79  FB4600 


80  N39825 


81  N63042 


82  N00204 


83  W22PL1 


Ft  Rucker 

AL 

Ft  Stewart 

GA 

Anniston 

AL 

Eglin  AFB 

FL 

Chicago 

IL 

San  Antonio 

TX 

Anniston 

AL 

Columbus 

OH 

Stockton 

CA 

Ft  Hood 

TX 

Monterrey 

CA 

San  Diego 

CA 

Ft  Benning 

GA 

Shaw  AFB 

SC 

Ft  Campbell 

KY 

Norfolk 

VA 

Philadelphia 

PA 

Tyndall  AFB 

FL 

Ft  Meade 

MD 

Jacksonville 

FL 

Ft  Huachuca 

AZ 

Tobyhanna 

PA 

Colts  Neck 

NJ 

Corpus  Christ! 

TX 

Alameda 

CA 

New  York 

NY 

Ft  Lewis 

WA 

Richmond 

VA 

Peterson  AFB 

CO 

North  Charleston 

SC 

Oakland 

CA 

Kessler  AFB 

MS 

Towson 

MD 

Oakland 

CA 

Nellis  AFB 

NV 

Huntsville 

AL 

Offutt  AFB 

NE 

Goose  Creek 

SC 

Lemoore  NAS 

CA 

Pensacola 

FL 

Ft  Knox 

KY 

Scott  AFB 

IL 

Philadelphia 

PA 

Oak  Harbor 

WA 

Fort  Worth 

TX 

Bayonne 

NJ 

Ft  Gillem 

GA 

Ft  Lee 

VA 

182,473 


180,196 


160,509 


158,219 


158,156 


149,595 


148,424 


144,464 


143,486 


141,391 


137,260 


135,608 


130,678 


130,099 


125,987 


117,907 


117,374 


115,670 


114,313 


114,279 


114,151 


113,712 


112,748 


108,233 


108,109 


107,829 


104,958 


102,938 


101,766 


99,517 


99,325 


99,316 


99,274 


97,741 


93,946 


91,313 


84,936 


84,201 


83,224 


81,383 


79,037 


78,149 


75,883 


75,506 


74,920 


72,379 


72,150 


236 


386 


129 


362 


99 


54 


159 


6,499 


806 


49 


91 

N00109 

Yorktown 

VA 

72,000 

51 

FB4626 

Malmstrom 

MT 

71,730 

68 

93 

SW3216 

Bremerton 

WA 

71,123 

386 

TOTAL  WEIGHT  (lbs.)  SHIPMENTS 

Top  93  39,959,720  163,803 

tons  19,980 

System  47,228,046  184,200 

tons  23,614 

PERCENT  84.61  88.93 

50 


Appendix  B:  Aerial  Ports  of  Debarkation 
(1  April  -  30  June  1997) 


APOD/ATIC 

CITY 

COUNTRY 

WEIGHT  (lbs.) 

1 

RMS 

RAMSTEIN 

GERMANY 

15,698,706 

2 

OSN 

OSAN 

KOREA  SOU 

3,681,797 

3 

HIK 

HONOLULU 

UNITED  ST 

2,183,901 

4 

OKO 

TOKYO 

JAPAN 

2,113,672 

5 

KWI 

KUWAIT 

KUWAIT 

1,945,244 

6 

HOW 

HOWARD 

PANAMA 

1,567,047 

7 

DHA 

DHAHRAN 

SAUDI  ARA 

1,452,686 

8 

DNA 

KADENA 

JAPAN 

1,261,717 

SIZ 

SIGONELLA 

ITALY 

1,180,224 

10 

MHZ 

MILDENHALL 

UNITED  Kl 

1,060,807 

KEF 

KEFLAVIK 

ICELAND 

1,033,173 

12 

BAH 

BAHRAIN 

BAHRAIN 

992,246 

13 

NBW 

GUANTANAMO 

CUBA 

737,502 

14 

NRR 

ROOSEVELT  ROADS 

PUERTO  RICO 

726,723 

15 

UAM 

ANDERSON  AFB 

GUAM 

629,812 

THU 

THULE  AIR  BASE 

GREENLAND 

579,544 

RTA 

ROTA  (NAS) 

SPAIN 

504,899 

18 

EDF 

ANCHORAGE 

UNITED  ST 

472,560 

19 

RUH 

RIYADH 

SAUDI  ARA 

432,963 

KWAJALEIN 

TRUST  TER 

398,640 

21 

PLA 

PALMEROLA 

HONDURAS 

389,241 

22 

ASP 

ALICE  SPRINGS 

AUSTRALIA 

378,266 

23 

PAP 

PORTAU  PRINCE 

HAITI 

274,966 

24 

TZL 

TUZLA 

BOSNIA-HE 

258,529 

MSJ 

MISAWA 

JAPAN 

243,966 

26 

EIL 

FAIRBANKS 

UNITED  ST 

224,657 

27 

LPB 

LAPAZ 

BOLIVIA 

215,073 

28 

LGS 

LAJES 

PORTUGAL 

191,941 

29 

ADA 

ADANA 

TURKEY 

190,497 

30 

UlO 

QUITO 

ECUADOR 

184,008 

31 

NKW 

DIEGO  GARCIA  (SEEA/0 

BR.  IND. 

164,186 

32 

TTH 

THUMRAIT 

OMAN 

150,224 

SAL 

SAN  SALVADOR 

EL  SALVAD 

144,976 

34 

JON 

JOHNSTON  ATOLL 

JOHNSTON 

142,051 

35 

FUK 

FUKUOKA 

JAPAN 

120,975 

TLV 

TEL  AVIV 

ISRAEL 

114,005 

ROM 

RICHMOND 

AUSTRALIA 

112,744 

ADH 

ALDAN 

RUSSIAN  F 

111,705 

MIQ 

CARACAS  (NOT  lATA;  S 

VENEZUELA 

105,112 

40 

SGP 

SINGAPORE 

SINGAPORE 

96,553 

BKK 

BANGKOK 

THAILAND 

90,748 

DKR 

DAKAR 

SENEGAL 

89,438 

SHIPMENTS 


36,528 

8,986 


10,408 


3,728 

474 

415 

131 

521 


659 

594 

91 

1,184 

2,567 

277 

236 


51 


43  NBO 


44  |KUZ 


46  IFUJ 


47  IWA 


48  pZP 


49  AVB 


51  UTP 


52 


53  ITIF 


56  TGU 


59  IKHE 


60  GUA 


61  RIO 


62  SIX 


64  IBUE 


65  |SCL 


67  CAI 


68  |SOC 


70  ICHC 


72  lAJR 


75  loco 


77  PSE 


79  NDJ 


80  FNA 


81  KWJ 


82  pJK 


83  ESB 


84  plY 


85  PBM 


86  SOX 


87  TJS 


88  BZE 


89  YOD 


90  IGL 


NAIROBI 

KENYA 

83,643 

KUNSAN 

KOREA  SOU 

82,960 

NAPLES 

Ftaly 

81,784 

FUKUE 

JAPAN 

77,707 

IWAKUNI 

JAPAN 

75,103 

SEVILLA 

SPAIN 

71,129 

AVIANO 

ITALY 

70,542 

WAKE  ISLAND 

WAKE  ISLA 

63,979 

RAYONG 

THAILAND 

59,378 

TASZAR 

HUNGARY 

58,894 

TAIF 

SAUDI  ARA 

53,483 

LIMA-CALLOA 

PERU 

51,480 

BOGOTA 

COLUMBIA 

50,833 

TEGUCIGALPA 

HONDURAS 

48,507 

AKROTIRI 

48,374 

WOOMERA 

AUSTRALIA 

46,025 

KIMHAE 

KOREA  SOU 

42,191 

GUATEMALA  CITY 

GUATEMALA 

33,681 

RIO  DE  JANIERO 

BRAZIL 

33,500 

ST  CROIX 

VIRGIN  IS 

31,602 

OLBIA 

ITALY 

29,895 

BUENOS  AIRES 

ARGENTINA 

29,824 

SANTIAGO 

CHILE 

29,079 

BRASILIA 

BRAZIL 

28,220 

CAIRO 

EGYPT 

28,162 

SOLO 

INDONESIA 

26,427 

AMMAN 

JORDAN 

26,340 

CHRISTCHURCH 

NEWZEALA 

25,742 

CUBI  POINT  BATAAN 

PHILIPPIN 

23,500 

ARVIDSJAUR 

SWEDEN 

22,900 

KAPIT 

MALAYSIA 

21,155 

MANAGUA 

NICARAGUA 

^  20,492 

SAN  JOSE  OCCIDENTAL 

PHILIPPIN 

19,074 

MONTEVIDEO 

URUGUAY 

18,500 

PONCE 

PUERTO  Rl 

18,430 

KINSHASA 

ZAIRE 

18,324 

N'DJAMENA 

CHAD 

18,068 

FREETOWN 

SIERRA  LE 

17,200 

KWANGJU 

KOREA  SOU 

16,565 

JAKARTA 

INDONESIA 

15,918 

ANKARA 

TURKEY 

13,535 

DIYARBAKIR 

TURKEY 

13,143 

PARAMARIBO 

SURINAME 

11,710 

SHAHEED  MWAFFAQ  AB 

JORDAN 

11,510 

TANJUNG  SELOR 

INDONESIA 

10,530 

BELIZE 

BELIZE 

9,880 

COLD  LAKE 

CANADA 

9,460 

[ZMIR 

TURKEY 

9,066 
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91  ADQ 


92  YHM 


93  NIM 


95  IKIN 


96  SGZ 


97  SYA 


98  GEO 


99  FRF 


102  SKJ 


103  YQX 


104  NAS 


105  BRI 


107  BRW 


109  BDS 


110  KUL 


114  KER 


115  MDY 


116  PPG 


117  IFR 


118  ASM 


119  ANC 


120  SDQ 


121  TRS 


122  KPO 


123  QUI 


124  BTI 


125  HRG 


KODIAK 

UNITED  ST 

6,427 

HAMILTON 

CANADA 

6,320 

NIAMEY 

NIGER 

5,856 

ISTRES 

FRANCE 

5,204 

KINGSTON 

JAMAICA 

4,245 

SONGKHLA 

THAILAND 

3,126 

SHEMYA 

UNITED  ST 

2,940 

GEORGETOWN 

GUYANA 

2,780 

FRANKFURT 

GERMANY 

2,763 

PATRAS 

GREECE 

2,382 

SAMARINDA 

INDONESIA 

2,099 

SITKINAK  ISLAND 

UNITED  ST 

2,000 

GANDER 

CANADA 

1,584 

NEW  PROVIDENCE  ISLAN 

BAHAMA  IS 

1,100 

BARI 

ITALY 

1,019 

SARAJEVO 

BOSNIA-HE 

900 

BARROW 

UNITED  ST 

633 

MARMANDE 

FRANCE 

600 

BRINDISI 

ITALY 

485 

KUALA  LUMPUR  INTL 

MALAYSIA 

391 

ISTANBUL 

TURKEY 

352 

DOHA 

QATAR 

208" 

ZAGREB 

CROATIA 

161 

KERMAN 

IRAN 

151 

MIDWAY  ISLAND 

MIDWAY  IS 

86 

PAGO  PAGO 

AMERICAN 

75 

ISTRES 

FRANCE 

67 

ASMARA 

ERITREA 

45 

ANCHORAC3E 

UNITED  ST 

26 

SANTO  DOMINGO 

DOMINICAN 

23 

TRIESTE 

ITALY 

17 

POHANG 

KOREA  SOU 

16 

CHUQUICAMATA 

CHILE 

7 

BARTER  ISLAND 

UNITED  ST 

6 

HURGHADA 

EGYPT 

5 

PALMA  DE  MALLORCA 

SPAIN 

3 

MUSCAT 

OMAN 

1 

TOTAL 

WEIGHT  (lbs.) 

SHIPMENTS 

top  22 

39,421,370 

160,579 

tons  19,711 

system 

44,409,271 

183,215 

tons  22,205 

PERCENT 

88.77 

87.65 
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Appendix  C:  Origin-APOD  Demand  Matrix. 
(Weight  in  Pounds  from  1  April  -  30  June  1996) 
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Appendix  D:  Macro  to  Produce  CPLEX  Readable  Linear  Program 

'  The  following  set  of  subroutines  creates  a  CPLEX  readable  file 
'  containing  the  transshipment  problem  examined  in  the  1997 
'  thesis  of  Capt  L  Dingle. 


Macro  programmed  by:  Maj  R  Hill,  AFIT/ENS 


Cautionary  notes: 

This  macro  is  not  intended  to  be  a  general  purpose  formulation 
tool.  As  such  there  are  certain  aspects  of  the  code  that  look  for 
specific  data  in  specific  locations  in  the  Main  worksheet. 
Furthermore,  input  error  checking  is  kept  to  a  minimum. 


Declare  global  /  public  variables 


Public  OriginNames ( 1  To  60)  As  String 
Public  APOENamesd  To  15)  As  String 
Public  APODNamesd  To  25)  As  String 
Public  Supply (1  To  60)  As  Double 
Public  Demand (1  To  25)  As  Double 
Public  APOECapacity (1  To  15)  As  Double 
Public  APOEFixedCosts { 1  To  25)  As  Double 
Public  APOEThruPutCosts (1  To  25)  As  Double 


Sub  OutputCPLEX { ) 
Sheets ( "Main” ) . Select 


Open  desired  output  file  in  which  to  place  the  data 


sPath  =  Cells (8,  15) 
sFile  =  Cells (10,  15) 

xTarget  =  sPath  &  "\"  &  sFile  &  ".dat" 
yTarget  =  sPath  &  "\"  &  sFile  &  ".idx" 
Open  xTarget  For  Output  As  #10 
windex  =  Cells (12,  15) 

If  windex  =  "yes"  Then 
IndexFlag  =  True 
End  If 


Next  read  the  Origin,  APOE,  and  APOD  data  into  the  public 
storage  arrays.  Keeping  the  data  in  memory  speeds  up  the 
processing  as  opposed  to  accessing  the  cells  directly  from  the 
spreadsheet . 


NumOf Origins  =  Cells (3,  1) 

NumOfAPOEs  -  Cells (16,  11) 

NumOfAPODs  =  Cells (3,  6) 

Read  in  the  origin  data,  names  and  supply  information 
This  data  starts  in  row  four  of  the  sheet  "Main" 
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For  i  =  1  To  NumOfOrigins 

Supply (i)  =  Cells (i  +  3,  3) 

OriginNames (i)  =  Cells (i  +  3,  2) 

Next  i 

TotalSupply  =  Cells (3,  3) 

Read  in  the  APOD  data,  names  and  demand  information 


For  i  =  1  To  NumOfAPODs 

Demand(i)  =  Cells (i  +  3,  9) 
APODNames(i)  =  Cells (i  +  3,  8) 
Next  i 

TotalDemand  =  Cells (3,  9) 


The  BigM  value  is  used  to  spoof  the  throughput  constraints  on 
each  of  the  APOEs.  Currently  there  is  no  constraint  on  throughput. 
However,  we  wish  to  have  that  capability  built  into  the  formulation. 

BigM  =  Application. Max (TotalSupply,  TotalDemand) 

Read  in  the  APOE  data,  names,  Capacities,  Operating  Costs 


For  i  =  1  To  NumOfAPOEs 

APOENames(i)  =  Cells (i  +  16,  11)  &  "  —  "  &  Cells (i  +  16,  13) 
APOEFixedCosts (i)  =  Cells (i  +  16,  15) 

APOEThruPutCosts (i)  =  Cells (i  +  16,  16) 

If  Cells (i  +  16,  14)  =  0  Then 
APOECapacity (i)  =  BigM 

Else 

APOECapacity (i)  =  Cells (i  +  16,  14) 

End  If 

Next  i 

Next  read  in  the  shipping  cost  data 

CostFromOrigin  =  Cells (2,  15) 

CostFromAPOE  =  Cells (4,  15) 


*  The  following  code  generates  an  index  of  the  Origin,  APOE,  and 

'  APOD  names  according  to  the  name  used  within  the  CPLEX  formulation. 

’  The  index  file  name  matches  the  formulation  name  with  the  exception 

'  of  using  the  .idx  suffix  versus  the  .dat  suffix. 

! 

If  IndexFlag  Then 

Open  yTarget  For  Output  As  #9 

Print  #9,  "Variables  in  this  model  are  of  the  following  form:" 
Print  #9,  "  0#A#  —  Flow  from  Origin  node  #  to  APOE  node 

#" 

Print  #9,  "  A#D#  —  Flow  from  APOE  node  #  to  APOD  node  #" 

Print  #9,  "where  the  node  numbers  are  defined  in  the  following 

manner : " 

Print  #9,  "First  is  the  list  of  Origin  nodes  in  the  model:" 

For  i  =  1  To  NumOfOrigins 

Print  #9,  "Origin  node  ",  i,  "  coded  as  ",  "0"  & 

LTrim(Str (i) ) ,  "  is  ",  OriginNames (i) 

Next  i 

Print  #9,  "Next  is  the  list  of  APOE  nodes  in  the  model:" 

For  i  =  1  To  NumOfAPOEs 
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coded  as 


»r  II  7\  II 


Print  #9,  "APOE  node  ”,  i,  " 

LTrim (Str (i) ) ,  ”  is  ",  APOENames (i ) 

Next  i 

Print  #9,  "Finally  is  the  list  of  APOD  nodes  in  the  model:" 
For  i  -  1  To  NumOfAPODs 

Print  #9,  "APOD  node  ",  i,  "  coded  as  ",  "D"  & 

LTrim (Str (i) ) ,  "  is  ",  APODNames(i) 

Next  i 

Print  #9,  "Variable  Z(J)  is  0  if  APOE  J  closed,  1  if  APOE  J 


IS  open" 
D#" 

labeled. " 


constraint" 
constraint " 


Print  #9,  "Variable  YA#D#  is  1  if  this  APOE  uniquely  serves 
Print  #9,  "  " 

Print  #9,  "Each  of  the  constraints  in  the  formulation  is 

Print  #9,  "The  following  is  the  coding  used  for  the  labels:" 

Print  #9,  "  Ob j :  -  Objective  function  of  the  problem" 

Print  #9,  "  S# :  -  Supply  constraints" 

Print  #9,  "  D# :  -  Demand  constraints" 

Print  #9,  "  FC# :  -  Net  flow,  or  flow  capacity 

Print  #9,  "  UJ# :  -  Unique  APOE  to  APOD  junction 

Print  #9,  "  PC#:  -  Aerial  port  capacity  constraint" 

Print  #9,  "  LC# :  -  Link  capacity  to  APOD  constraints" 

Print  #9,  "  SC#:  -  Special  demand  constraints" 

Print  #9,  "End  of  List...." 


Close  #9 
End  If 


*  The  following  section  of  the  macro  formats  and  outputs  the 
objective 

*  function  for  the  problem. 


'  There  are  three  pieces  of  the  objective  function.  The  first 

piece 

^  represents  the  cost  of  shipping  material  from  the  Origin  to  the 
APOE. 

*  The  second  piece  represents  the  cost  of  shipping  material  from  the 
APOE 

*  to  the  APOD.  The  final  piece  represents  the  fixed  operating  costs 
for 

*  opening  up  an  APOE . 


Print  #10,  "Minimize" 


TermsPerLine  =  5 

Sheets { "OriginDistances" ) . Select 
xRow  =  4  +  NumOfOrigins 
yCol  =  2  +  NumOfAPOEs 

I 

*  Create  the  first  piece  of  the  objective  function. .. the  Origin  to 
APOE 

'  shipment  costs.  The  code  includes  test  for  zero  objective  function 
values , 

*  which  are  skipped  over. 

» 

Set  xRange  =  ActiveSheet . Range (Cells ( 5,  3),  Cells (xRow,  yCol) ) 
WorkTerm  =  "Obj :  " 
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Numlnit  -  0 
For  Each  C  In  xRange 
xRow  =  C . Row 
yCol  =  C. Column 

zCost  =  CostFromOrigin  *  C. Value  +  APOEThruPutCosts (Cells (3, 

yCoD) 

If  zCost  <>  0  Then 

WorkTerm  =  WorkTerm  &  LTrim (Str  (zCost )  )  &  "  ** 

WorkTerm  =  WorkTerm  &  "0"  &  LTrim (Str (Cells (xRow,  1))) 
WorkTerm  =  WorkTerm  &  "A"  &  LTrim (Str (Cells ( 3 ,  yCol) ) )  &  "  + 

IT 

Numlnit  =  Numlnit  +  1 
End  If 

If  Numlnit  >=  TermsPerLine  Then 
Print  #10,  WorkTerm 
WorkTerm  = 

Numlnit  =  0 
End  If 
Next  C 


'  Create  the  second  piece  of  the  objective  function. .. the  APOE  to 

APOD 

'  shipment  costs. 

*  Notice  I  do  not  want  to  reset  the  WorkTerm  or  Numlnit  since  I  am 
still 

*  building  the  same  objective  function  value. 

^  The  code  includes  test  for  zero  objective  function  values, 

*  which  are  skipped  over. 


Sheets ( "APODDistances ” ) .Select 
xRow  =  4  +  NumOfAPODs 
yCol  =  2  +  NumOfAPOEs 

Set  xRange  =  ActiveSheet . Range (Cells (5,  3),  Cells (xRow,  yCol) ) 

For  Each  C  In  xRange 
xRow  =  C.Row 
yCol  “  C. Column 

zCost  =  CostFromAPOE  *  C. Value 
If  zCost  <>  0  Then 

WorkTerm  =  WorkTerm  &  LTrim (Str ( zCost ) )  &  ”  ” 

WorkTerm  =  WorkTerm  &  "A"  &  LTrim (Str (Cells (2,  yCol) ) ) 
WorkTerm  =  WorkTerm  &  ”D"  &  LTrim (Str (Cells (xRow,  1)))  &  "  + 

It 

Numlnit  =  Numlnit  +  1 
End  If 

If  Numlnit  >=  TermsPerLine  Then 
Print  #10,  WorkTerm 
WorkTerm  = 

Numlnit  =  0 
End  If 
Next  C 


*  Create  the  last  piece  of  the  objective  function ...  the  APOE  fixed 
’  operating  costs  given  that  the  APOE  has  been  opened.  The  decision 
to  open 

’  the  APOE  is  captured  in  the  Z#  variable  within  the  formulation. 

I 

'  Notice  I  still  do  not  want  to  reset  the  WorkTerm  or  Numlnit 

since  I 

*  building  the  same  objective  function  value. 
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For  j  =  1  To  NumOfAPOEs  “  1 

WorkTerm  =  WorkTerin  &  LTrim  (Str  (APOEFixedCosts  ( j  )  )  )  &  "  Z"  & 
LTrim(Str { j ) )  &  ”  +  " 

Numlnlt  =  Numlnit  +  1 
If  Numlnlt  >=  TermsPerLine  Then 
Print  #10,  WorkTerm 
WorkTerm  =  "" 

Numlnlt  =  0 
End  If 
Next  j 

WorkTerm  =  WorkTerm  &  LTrim (Str (APOEFixedCosts (NumOfAPOEs )) ) 
WorkTerm  =  WorkTerm  &  "  Z”  &  LTrim (Str (j ) ) 

Print  #10,  WorkTerm 

I 

*  Now  ready  to  start  formating  and  printing  out  the  constraints 

t 

’  All  but  the  Special  Demand  constraints  are  generated  based  on 

the  data 

*  stored  in  the  public  arrays  previously  filled  from  the  "MAIN" 
worksheet . 

t 

TermsPerLine  =  8 
Print  #10,  "Subject  to" 

I 

*  The  next  section  of  code  will  generate  the  Supply  constraints 

*  Currently  the  code  assumes  that  each  origin  will  have  non-zero 
supply. 

*  Thus  there  is  no  checking  for  zero  values. 

I 

’Print  #10,  "Following  are  the  Supply  constraints:" 

ConstraintNumber  =  1 
For  i  =  1  To  NumOfOrigins 

WorkTerm  =  "S"  &  LTrim (Str (ConstraintNumber ) )  &  ":  " 
ConstraintNumber  =  ConstraintNumber  +  1 
For  j  =  1  To  NumOfAPOEs  -  1 

WorkTerm  =  WorkTerm  &  "0"  &  LTrim (Str (i ) )  &  "A"  &  LTrim (Str (j ) ) 

&  "  +  ” 

Next  j 

WorkTerm  =  WorkTerm  &  "0"  &  LTrim (Str (i ) )  &  "A"  & 

LTrim (Str (NumOfAPOEs) )  &  "  <=  " 

WorkTerm  =  WorkTerm  &  Str ( Supply (i ) ) 

Print  #10,  WorkTerm 
Next  i 


’  The  next  section  of  code  will  generate  the  Demand  constraints 

’  As  with  the  supply  constraints,  the  current  code  that  follows 
assumes 

’  that  each  APOD  will  have  a  non-zero  demand. 

t 

'Print  #10,  "Following  are  the  Demand  constraints:" 

ConstraintNumber  =  1 
For  i  =  1  To  NumOfAPODs 

WorkTerm  =  "D"  &  LTrim (Str (ConstraintNumber ) )  &  ":  " 
ConstraintNumber  =  ConstraintNumber  +  1 
For  j  =  1  To  NumOfAPOEs  -  1 

WorkTerm  =  WorkTerm  &  "A"  &  LTrim (Str (j ) )  &  "D"  &  LTrim (Str (i) ) 

&  "  +  " 

Next  j 

WorkTerm  =  WorkTerm  &  "A"  &  LTrim (Str (NumOfAPOEs ) )  &  "D"  & 

LTrim (Str (i) )  &  »  =  " 

WorkTerm  =  WorkTerm  &  Str (Demand (i) ) 
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Print  #10,  WorkTerm 
Next  i 


'  The  next  section  of  code  will  generate  the  Flow  Conservation 

constraints . 

’  Essentially  a  flow  conservation  constraint  ensures  that  all  goods 
that  flow 

’  into  an  APOE  in  fact  flow  out ...  In-flow  =  Out-flow. 


*  Print  #10,  "Following  are  the  Flow  Conservation  constraints:" 
ConstraintNumber  =  1 
For  j  =  1  To  NumOfAPOEs 
Termsinit  =  0 

WorkTerm  =  "FC"  &  LTrim (Str (ConstraintNumber) )  &  " 

ConstraintNumber  =  ConstraintNumber  +  1 
For  i  =  1  To  NumOfOrigins  -  1 

WorkTerm  =  WorkTerm  &  "0"  &  LTrim (Str (i ) ) 

WorkTerm  =  WorkTerm  &  "A"  &  LTrim (Str (j ) )  &  "  +  " 
Termsinit  =  Termsinit  +  1 
If  Termsinit  >=  TermsPerLine  Then 
Print  #10,  WorkTerm 
WorkTerm  =  "" 

Termsinit  =  0 


End  If 
Next  i 

WorkTerm  =  WorkTerm  &  "0"  &  LTrim (Str (NumOfOrigins ) ) 
WorkTerm  =  WorkTerm  &  "A"  &  LTrim (Str (j ) )  &  "  -  " 
Termsinit  ==  Termsinit  +  1 
For  k  =  1  To  NumOfAPODs  -  1 

If  Termsinit  >=  TermsPerLine  Then 
Print  #10,  WorkTerm 
WorkTerm  =  "" 

Termsinit  =  0 
End  If 


WorkTerm  =  WorkTerm  & 
WorkTerm  =  WorkTerm  & 
Termsinit  =  Termsinit 
Next  k 

WorkTerm  =  WorkTerm  &  "A" 
WorkTerm  =  WorkTerm  &  "D" 
Print  #10,  WorkTerm 
Next  j 


"A"  &  LTrim ( Str (j) ) 

"D"  &  LTrim { Str (k) )  &  "  -  " 

+  1 

&  LTrim(Str (j ) ) 

&  LTrim (Str (NumOfAPODs) )  &  " 


0.0" 


The  next  section  of  code  will  generate  the  constraints  ensuring 
that  each  APOE  uniquely  serves  an  APOD. 


’Print  #10,  "Following  are  the  APOE-APOD  uniqueness  constraints:" 
ConstraintNumber  =  1 
For  k  =  1  To  NumOfAPODs 

WorkTerm  =  "UJ"  &  LTrim (Str (ConstraintNumber ) )  &  ":  " 
ConstraintNumber  =  ConstraintNumber  +  1 
For  j  =  1  To  NumOfAPOEs  -  1 

WorkTerm  =  WorkTerm  &  "YA"  &  LTrim (Str (j ) )  &  "D"  &  LTrim (Str ( k) ) 

&  "  +  " 

Next  j 

WorkTerm  =  WorkTerm  &  "YA"  &  LTrim (Str (NumOfAPOEs) )  &  "D"  & 
LTrim(Str (k) )  &  "  =  1" 

Print  #10,  WorkTerm 
Next  k 
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*  The  next  section  of  code  will  generate  the  Aerial  port  capacity 

’  constraints.  Currently,  there  is  no  requirement  to  enter  capacities 

*  on  any  of  the  APOEs.  When  no  capacity  is  specified,  this  macro  uses 

*  the  BigM  value  (method)  for  ensuring  sufficient  capacity  for  each 
APOE. 

*  BigM  is  taken  here  as  the  maximum  of  Total  Supply  or  Total 
Demand. 


’Print  #10,  "Following  are  the  Aerial  Port  capacity  constraints:" 
ConstraintNumber  =  1 
^  For  j  =  1  To  NumOfAPOEs 

WorkTerm  ==  "PC"  &  LTrim  (Str  (ConstraintNumber )  )  &  " 

ConstraintNumber  =  ConstraintNumber  +  1 
Termsinit  =  0 

For  i  =  1  To  NumOfOrigins  -  1 

WorkTerm  =  WorkTerm  &  "0"  &  LTrim (Str (i) )  &  "A"  &  LTrim (Str (j ) ) 

&  "  +  " 

Termsinit  =  Termsinit  +  1 
If  Termsinit  >=  TermsPerLine  Then 
Print  #10,  WorkTerm 
WorkTerm  =  "" 

Termsinit  =  0 
End  If 
Next  i 

WorkTerm  =  WorkTerm  &  "0"  &  LTrim (Str (NumOfOrigins ) )  &  "A"  & 

LTrim (Str ( j ) ) 

WorkTerm  =  WorkTerm  &  "  -  "  &  LTrim (Str (APOECapacity (j )) )  &  "  Z"  & 

LTrim ( Str { j ) ) 

WorkTerm  =  WorkTerm  &  "  <=  0.0  " 

Print  #10,  WorkTerm 
Next  j 


'  The  following  section  generates  APOE  to  APOD  link  capacity 

constraints . 

'  These  constraints  kick  in  when  a  given  APOE  to  APOD  link  is 
initialized  to 

'  uniquely  serve  each  of  the  APOD. 

’  There  is  one  constraint  for  each  of  the  APOE  ~  APOD 

combinations . 

r 

’  Here  again,  there  is  no  specification  of  link  capacities  within 

the 

’  provided  data.  For  this  coding,  the  logic  is  to  ensure  that  no  link 

’  ships  more  than  the  capacity  of  the  APOE  from  which  it  originates. 

Thus 

’  the  BigM  value  is  again  used  in  lieu  of  true  capacity  constraint 
value . 


’Print  #10,  "Following  are  the  Link  capacity,  APOE  to  APOD  constraints:" 
ConstraintNumber  =  1 
For  j  =  1  To  NumOfAPOEs 

For  k  =  1  To  NumOfAPODs 

WorkTerm  =  "LC"  &  LTrim (Str (ConstraintNumber ) )  &  ":  " 

ConstraintNumber  =  ConstraintNumber  +  1 

WorkTerm  =  WorkTerm  &  "A"  &  LTrim (Str (j ) )  &  "D"  &  LTrim (Str ( k) ) 
WorkTerm  =  WorkTerm  &  "  -  "  &  LTrim (Str (APOECapacity (j )) )  &  "  " 
WorkTerm  =  WorkTerm  &  "YA"  &  LTrim (Str (j ) )  &  "D"  &  LTrim ( Str ( k) ) 
WorkTerm  =  WorkTerm  &  "  <=  0.0" 
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Print  #10,  WorkTerm 
Next  k 
Next  j 


The  following  section  generates  the  special  consideration  supply 
to 

*  demand  constraints  based  on  the  data  within  the  Special  Demand 
matrix 

*  contained  in  the  SpecialDemand  sheet. 

I 

1 

’  Define  the  special  demand  matrix  as  an  object  and  then  examine 

*  each  cell  in  the  defined  range. 

t 

*  Print  #10,  "Following  are  the  Special  Demand  constraints:" 

Sheets ( "SpecialDemand" ) . Select 

xRow  =  6  +  NumOfOrigins 
yCol  =  3  +  NumOfAPODs 

Set  xRange  =  ActiveSheet . Range (Cells (7,  4),  Cells (xRow,  yCol) ) 

'  Next  consider  each  cell  in  the  defined  range 


ConstraintNumber  =  1 
For  Each  C  In  xRange 

If  C. Value  >  0  Then 
xRow  =  C.Row, 
yCol  =  C. Column 
For  j  =  1  To  NumOfAPOEs 

WorkTerm  =  "SC"  &  LTrim (Str (ConstraintNumber ) ) 
ConstraintNumber  =  ConstraintNumber  +  1 


WorkTerm  =  WorkTerm  &  "0"  & 
WorkTerm  =  WorkTerm  &  "A"  & 
WorkTerm  =  WorkTerm  &  LTrim 
WorkTerm  =  WorkTerm  &  "  YA" 
WorkTerm  =  WorkTerm  &  "D"  & 
WorkTerm  =  WorkTerm  &  "  >=  0 
Print  #10,  WorkTerm 
Next  j 
End  If 
Next  C 


LTrim (Str (Cells (xRow,  1) ) ) 
LTrim(Str (j) )  &  "  -  " 

Str (C. Value) ) 

&  LTrim (Str ( j ) ) 

LTrim(Str (Cells (4,  yCol) ) ) 
0" 


*  All  the  standard  constraints  are  now  placed  into  the  file. 

'  The  following  section  adds  the  non-negativity  constraints  and 

'  the  binary  constraints  to  the  problem. 

? 

'  The  default  bounds  are  0  <=  x  <=  infinity.  These  default  bounds 

'  apply  to  each  of  the  0#A#  and  A#D#  variables  in  the  problem. 

I 

*  In  this  initial  formulation,  the  YA#D#  variables  are  allowed  to 
be 

’  real  variables  bounded  above  by  1 .  In  actuality,  these  will  likely 
have 

'  to  be  recoded  as  Integer  variables. 


Print  #10,  "Bounds" 

For  j  =  1  To  NumOfAPOEs 

For  k  =  1  To  NumOfAPODs 

WorkTerm  =  "YA"  &  LTrim (Str (j ) )  &  "D"  &  LTrim (Str (k) ) 
WorkTerm  =  WorkTerm  &  "  <=  1" 
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Print  #10,  WorkTerm 
Next  k 
Next  j 


*  Integer  section  follows.  By  default,  integer  variables  are 

assumed 

’  binary  and  provided  the  bounds  of  0<=  y  <=  1.  This  assumption  is 
used. 

f 

Print  #10,  "Integers” 

WorkTerm  =  "" 

For  j  =  1  To  NumOfAPOEs 

WorkTerm  =  WorkTerm  &  "Z"  &  LTrim (Str { j ) )  &  "  " 

Next  j 

Print  #10,  WorkTerm 


*  Finally  END  the  file 

I 

Print  #10,  "End" 

t 

*  Close  the  file  and  return  to  the  "MAIN"  worksheet 

I 

Close  #10 

Sheets ( "Main" ) .Select 

Beep 

Beep 

End  Sub 


Sub  Oops ( ) 

Close 
End  Sub 

Sub  auto_open ( ) 

Sheets ( "Modulel" ). Visible  =  False 
End  Sub 

Sub  Reshow  0 

Sheets  (  "Modulel" ).  Visible  ==  True 
End  Sub 
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